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Abstract 



We develop an operator formalism for investigating the properties of nonabelian 
• cosmic strings (and vortices) in quantum field theory. Operators are constructed 

that introduce classical string sources and that create dynamical string loops. The 
^ . operator construction in lattice gauge theory is explicitly described, and correlation 

functions are computed in the strong-coupling and weak-coupling limits. These cor- 
relation functions are used to study the long-range interactions of nonabelian strings, 
taking account of charge-screening effects due to virtual particles. Among the phe- 
nomena investigated are the Aharonov-Bohm interactions of strings with charged par- 
ticles, holonomy interactions between string loops, string entanglement, the transfer 
of "Cheshire charge" to a string loop, and domain wall decay via spontaneous string 
nucleation. We also analyze the Aharonov-Bohm interactions of magnetic monopoles 
with electric flux tubes in a confining gauge theory. We propose that the Aharonov- 
Bohm effect can be invoked to distinguish among various phases of a nonabelian gauge 
theory coupled to matter. 



1. Introduction 



In many cases, a gauge theory that undergoes the Higgs mechanism will con- 
tain topologically stable defects. Among the defects that can occur are line defects 
("cosmic strings") in 3+1 dimensions, and pointhke vortices in 2+1 dimensions!^' 

If the symmetry breaking pattern and the matter content of the theory are suit- 
able, the strings or vortices may have long range interactions with various particles. 
These interactions are a consequence of the Aharonov-Bohm phenomenon — the wave 
function of a particle acquires a nontrivial phase when the particle is covariantly 
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transported around the string. ' 

This simple observation has interesting implications!* Because the Aharonov- 
Bohm interaction has infinite range, and no local operator can destroy an object 
that has an infinite-range interaction with another object, gauge theories with such 
interactions always respect nontrivial superselection rules. The structure of the su- 
perselection sectors can be invoked to distinguish among the various possible phases 
of a gauge theory. Furthermore, the existence of infinite range interactions that are 
fundamentally quantum mechanical exposes the limitations of the "no-hair" theo- 
rems of black hole physics. Though a black hole may have no classical hair, it can 
carry quantum numbers that are detectable by means of the interaction between the 
black hole and a cosmic string!*''^ ^' 

The "Aharonov-Bohm phase" acquired by an object that circumnavigates a string 
is, in general, a gauge transformation contained in the unbroken gauge group. When 
the manifest gauge symmetry is nonabelian, this gauge transformation might not lie 
in the center of the group. In that case, we say that the string is "nonabelian." 
The physical properties of a nonabelian string are qualitatively different from the 
properties of an abelian string!^'*''" In particular, a loop of nonabelian string (or a 
pair of nonabelian vortices, in 2+1 dimensions) can carry a nontrivial gauge charge, 
so that the loop has an Aharonov-Bohm interaction with other strings. Moreover, 
nonabelian Aharonov-Bohm interactions involve transfer of charge between string 
loops and charged particles. Remarkably, the charge carried by a loop (which has a 
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topological origin) can not be localized anywhere on the string loop or in its vicinity. 
Following Ref. 13, we will refer to such unlocalized charge as "Cheshire charge." 

If we are to appeal to the Aharonov-Bohm effect to probe the phase structure of 
a gauge theory, or to investigate the quantum physics of black holes, we must be able 
to discuss interference phenomena in a language that does not rely on weak-coupling 
perturbative methods; we need a framework that (at least in principle) takes full 
account of the effects of virtual pairs and of the fluctuations of quantum fields. Such 
a framework was erected, for abelian strings (or vortices), in Ref. 5. There, operators 
were constructed that create a loop of string, or that introduce (as a classical source) 
the closed world sheet of a string. The correlation functions of these operators can 
be studied to investigate the properties of the strings, and their Aharonov-Bohm 
interactions in particular. 

Our main objective in this paper is to generalize the work of Ref. 5 to the case of 
nonabelian strings. Because of the subtle and elusive physics of nonabelian strings, 
this generalization is not entirely straightforward!"' Our primary motivation comes 
from two considerations. First, we seek assurance that the exotic physics of the 
nonabelian Aharonov-Bohm effect, previously inferred in the weak-coupling limit, 
actually survives in a fully quantum field-theoretic treatment. Second, we hope to 
construct (nonlocal) order parameters that can be used to classify the phases of a 
gauge theory. 

Let us formulate the order-parameter problem more explicitly, and in so doing, 
review some of the principal results of Ref. 5. 

A gauge theory can have an interesting phase diagram. Depending on its Higgs 
structure and on the parameters of the Higgs potential, the theory may be in a 
Coulomb (massless) phase, a Higgs phase, or a confinement phase. A Higgs phase is, 
roughly speaking, characterized by the existence of stable magnetic flux tubes, and a 
confinement phase by the existence of stable electric flux tubes. 

Nonlocal gauge-invariant order parameters can be devised that distinguish among 
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the various phases. The expectation value of the Wilson loop operator exhibits 
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area-law behavior if there are stable electric flux tubes, and perimeter-law behav- 
ior otherwise. The expectation value of the 't Hooft loop operator'^'"' exhibits area 
law behavior if there are stable magnetic flux tubes, and perimeter-law behavior 
otherwise. 

These order parameters are not sufficient, however, to distinguish among all pos- 
sible phases of a general gauge theory. Consider, for example, the case of an SU (N) 
gauge theory with matter in the fundamental representation. In this theory, the Wil- 
son loop always obeys the perimeter law, because an electric flux tube can break via 
nucleation of a pair. An 't Hooft loop operator can also be defined, but always obeys 
the area law.* Yet the theory can have a nontrivial phase diagram. Using adjoint 
Higgs fields, it is possible to break the gauge group down to its center Zjsf. This 
phase, which admits free Zn charges, is distinguishable from the confinement phase. 

In Ref. 5, an order parameter was described that can distinguish the free-charge 
phase from the confinement phase in an SU{N) gauge theory^ The idea is that the 
free-charge phase supports stable magnetic flux tubes (cosmic strings), and these 
strings have an infinite-range Aharonov-Bohm interaction with charges. No such 
interaction can exist if Zjv charges are conflned, or if Zn is spontaneously broken 
by a Higgs fleld that transforms as the fundamental representation. (Indeed, the 
confinement phase and the "Higgs" phase with Zn completely broken are not distin- 
guishable; this is an instance of "complementarity." 

To construct the order parameter, we must first devise an operator -F(S) that 
introduces a cosmic string world sheet on the closed two-dimensional surface S. (As 
we will discuss later, this operator is closely related to the 't Hooft loop operator.) 
Then consider the quantity 

(F(S))(iy(C)) ' ^^-^^ 



★ This is actually an oversimplification of the status of the 't Hooft loop, as we will discuss in 
Sections 5, 7, and 9. 

t This order parameter was actually discussed earlier by Predenhagen and Marcu in Ref. 17. See 
also Ref. 18. 
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where W{C) is the (fundamental representation) Wilson loop. In the free-charge 

phase, we have 



Here the limit is taken with E and C increasing to infinite size, and with the closest 
approach of E to C also approaching infinity; k^E, C) denotes the linking number of 
the surface E and the loop C. On the other hand, if there are no free Zjv charges, 
then we have 



The nonanalytic behavior of A(E, C) guarantees that a well-defined phase boundary 
separates the two phases. We will refer to A(Ti,C) as the "ABOP", or "Aharonov- 
Bohm Order Parameter." 

One of our objectives in this paper is to generalize the above construction, and 
to further explore its consequences. In general, it is not sufficient to determine the 
realization of the center of the gauge group, in order to distinguish all possible phases 
of a gauge theory. For example, SU (N) might break to a discrete subgroup H that is 
not contained in the center. Different unbroken groups (including nonabelian ones) 
can be distinguished according to the varieties of cosmic stings that exist, and the 
nature of the Aharonov-Bohm interactions of these strings with free charges. To 
probe the phase structure more thoroughly, we need to construct a more general F(E) 
operator. And we will need to consider carefully the interpretation of the behavior of 
the corresponding A(S, C) operator. The interpretation involves subtleties associated 
with the nonabelian Aharonov-Bohm effect. 

In the example described above, there are other order parameters that can serve 
to distinguish the free-charge phase from the confinement phase!"' But we believe 
that the properties of strings and of their Aharonov-Bohm interactions provide a more 
powerful method for classifying phases in more general cases. 

The remainder of this paper is organized as follows: Section 2 reviews the for- 
malism for describing configurations of many nonabelian vortices or strings. We 




(1.2) 



lim (A(S,C)) = 1 . 



(1.3) 



4 



emphasize the need to select an (arbitrary) basepoint for the purpose of defining 
the "magnetic flux" of a vortex, discuss the holonomy interactions between vortices, 
explain the origin of Cheshire charge, and note that nonabelian strings generically 
become entangled when they cross each other. 

In Section 3, we discuss the properties of domain walls that are bounded by loops 
of nonabelian string, and observe that physically distinct strings can bound physically 
identical walls. Hence when a wall "decays" by nucleating a loop of string, several 
competing "decay channels" may be available. 

Section 4 concerns Aharonov-Bohm interactions in abelian gauge-Higgs systems. 
We explore the screening of such interactions due to Higgs condensation. The obser- 
vations in Section 3 concerning domain wall decay are invoked, in order to interpret 
the results. 

In Section 5, we begin our analysis of the effects of quantum fluctuations on 
the nonabelian Aharonov-Bohm effect. We find that nonperturbative fluctuations 
cause ambiguities in the Aharonov-Bohm "phase" beyond those that occur in lead- 
ing semiclassical theory. We construct operators that create configurations of many 
nonabelian strings. Two types of operators are considered. One type introduces 
classical string sources on closed world sheets. The other type creates (or destroys) 
dynamical string loops. 

A subtle aspect of the construction is that, for both types, the many string con- 
figurations are coherent. This means the following: The "magnetic fiux" carried by 
a string can be characterized by an element a of the unbroken gauge group H. If 
an object that transforms as the representation (u) of H is transported about this 
string, the Aharonov-Bohm phase that it acquires, averaged over a basis for the 
representation, is 

-X^^ha) , (1.4) 

where n^, is the dimension of (u) and x^'^^ is its character!^'^^' Now if an a string and 
a b string are combined incoherently, then the averaged phase acquired by an object 
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that traverses the two strings in succession is 

-X^''\a) -x^''\b) ■ (1.5) 

that is, it is the product of the Aharonov-Bohm factors associated with the two 
individual strings. But if the two strings are patched together coherently, then the 
averaged phase becomes 

-X^^\ah) . (1.6) 

Til/ 

This property, that the Aharonov-Bohm factor associated with a pair of coherently 
combined strings is noi just the product of the Aharonov-Bohm factors of the two in- 
dividual strings, is a distinguishing feature of the nonabelian Aharonov-Bohm effect, 
and one of the goals of our work has been to better understand how this coherence is 
maintained when quantum fluctuations of fields are fully taken into account. 

Correlation functions of the operators constructed in Section 5 are used in Section 
6 to investigate Aharonov-Bohm interactions in a pure gauge theory, at both strong 
and weak coupling. In Section 7, we use these operators to study the quantum- 
mechanical mixing between different string states first discussed in Ref. 21. 

Holonomy interactions between pairs of vortices, or pairs of string loops, are con- 
sidered in Section 8. Correlation functions are used to demonstrate holonomy scatter- 
ing and string entanglement. We also analyze a correlation function that realizes the 

transfer of charge from a charged source to a loop of string, and the subsequent detec- 
tion of the (Cheshire) charge on the loop by means of its Aharonov-Bohm interaction 
with yet another loop of string. 

In Section 9, we consider nonabelian gauge-Higgs systems; we analyze the con- 
sequences of the Higgs mechanism for the stability of strings, and for the screening 
of Aharonov-Bohm interactions. We discuss how nonlocal order-parameters can be 
used to explore the phase structure of such a system. 

Section 10 concerns the properties of dynamical magnetic monopoles, in a confin- 
ing gauge theory. We construct generalizations of the Wilson and 't Hooft operators 
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for a theory with dynamical monopoles, and use these operators to demonstrate the 

[22 23] 

Aharonov-Bohm interaction between monopoles and electric flux tubes. 

The Appendix provides some additional details concerning some of the lattice 
calculations that are mentioned in the text of the paper. 



2. Nonabelian vortices and strings 

We will briefly review the formalism'^*' for describing conflgurations of many 
vortices (in two spatial dimensions) or many strings (in three spatial dimensions). 
Our purpose is to remind the reader of two peculiar features of nonabelian strings. 
First, there is an ambiguity when two or more loops of string are patched together 
to construct a multi-string configuration. Second, a loop of string can carry charges, 
and can exchange charge with other objects by means of the nonabelian Aharonov- 
Bohm effect. An understanding of these features is needed in order to interpret the 
behavior of the correlation functions that we will construct. 

2.1. Many- Vortex Configuration 

We will consider vortices first, and discuss strings later. A single isolated vortex 
can be associated with an element of the unbroken gauge group H, according to 







\ 


a(C, xq) = Pexp 




dx 




\ C,xo 


J 



Here C is a closed oriented path, far from the vortex core, that encloses the vortex and 
begins and ends at the point xq. This group element a(C, xq) is invariant under de- 
formations of the path C that keep xq fixed and that avoid the vortex core. An object 
that transforms as the irreducible representation (z/) of H acquires the "Aharonov- 
Bohm" phase D^'^\a{C,XQ)) when covariantly transported around the vortex. We 
require a & H, because the Higgs fields that drive the symmetry breakdown must 
return to their original values when so transported. The vortices can be topologically 
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classified, with the topological charge taking values in no{H); that is, a vortex with 
"flux" a E H can be smoothly deformed to another vortex with flux b E H if and 
only if a and b lie in the same connected component of H. If if is a discrete group 
(as we will usually assume in this paper), then, the topological charge is specified by 
an element of H. 

Similarly, we may associate n elements of H with a configuration of n vortices. 
To do so, we must choose n standard paths, all beginning and ending at the same 
point xq, that circumnavigate the various vortices!^"' 

This description of the n- vortex configuration evidently suffers from some ambi- 
guities. First, it is not gauge invariant. Under a gauge transformation that takes the 
value /i e if at the point xq (and so preserves the Higgs order parameter at that 
point), the elements ai, 02, . . . , On transform according to 

Qi haih~^ . (2.2) 

Thus, the configurations of a single vortex are labeled by the conjugacy classes of 
H. But the gauge freedom for many vortices involves just one overall conjugation. 
This means that specifying the positions and classes of all vortices is not sufficient to 
characterize the multi-vortex configuration. Single vortices can be patched together 
in various ways to construct different multi-vortex states. 

Another ambiguity concerns the choice of the paths Cj that enclose the vortices. 
The paths that begin and end at xq and avoid the cores of n vortices fall into homotopy 
classes; these classes are the elements of 7ri(M^,a;o), the fundamental group of Mn, 
the plane with n punctures. (This group is F^, the free group with n generators.) 
By assigning elements of H to each of the generators of 7ri(M„, xq), we have defined 
a homomorphism from 7ri(Mfj,xo) into H. This homomorphism (modulo the one 
overall conjugation) characterizes the configuration of n vortices. 

★ Here, we have implicitly defined H as the unbroken subgroup of a simply connected underlying 
(spontaneously broken) gauge group. The global topology of the underlying group is irrelevant, 
however, for the purpose of classifying the Aharonov-Bohm interactions of the vortices. 
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But this homomorphism is still ambiguous, because the n generators can be chosen 
in various ways. Consider in particular the case of two vortices, as shown inFig. 1. 
Standard paths ai and 0:2 have been chosen in Fig. la that wind counterclockwise 
around the two vortices. A topologically distinct choice for the path around vortex 
1 is shown in Fig. Ic, and a distinct choice for the path around vortex 2 is shown in 
Fig. Id. Suppose now that vortex 1 winds around vortex 2 (in the sense defined by 
the path a2), and returns to its original position, as in Fig. lb. We may deform our 
paths during the winding so that no vortex ever crosses any path; then each path is 
mapped to the same group element after the winding as before the winding. But after 
the winding, the final deformed path is not homotopically equivalent to the initial 
path. Therefore, the homomorphism that maps 7ri(M2, xq) to the group has changed. 

Suppose, for example, that initially ai is mapped to oi and a2 is mapped to 
a2- We wish to determine the final values, after the winding, of the group elements 
associated with the paths ai and a2- For this purpose, it is convenient to notice that, 
during the winding, the path shown in Fig. Ic is "dragged" to ai. Therefore, the 
group element associated with this path before the winding is the same as the group 
element associated with ai after the winding. We see that this path can be expressed 
as {aia2)ai{aia2)^^ , where ai«2 denotes the path that is obtained by traversing 
a2 first and ai second. The initial homomophism maps this path to a path-ordered 
exponential with the value 

a'l = (0102)01(0102)""^ . (2.3) 

Similarly we note that the path shown in Fig. Id is dragged during the winding to 
a2- This path is (q;iq;2)q;2(q!iQ!2)~^, and is mapped by the initial homomorphism to 

03 = (0102)02(0102)^-^ . (2.4) 

We conclude that the final homomorphism after the winding maps ai to a'^ and maps 
a2 to 02. Winding vortex 1 around vortex 2 has changed the two- vortex state (if 01 
and 02 do not commute) — both group elements have become conjugated by 01O2. 
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This change in the two-vortex state is not a mere mathematical curiosity. 
The physical interpretation is that there is a long-range interaction between non- 
commuting vorticcs!^ ''^^' Wc refer to this as the "holonomy interaction." It is a type of 
Aharonov-Bohm interaction, but we prefer not to use that name in this context. The 
holonomy interaction is really a classical effect, while the term "Aharonov-Bohm" is 
best reserved for an intrinsically quantum-mechanical interaction. When the number 
of vortices is three or more, these holonomy interactions can be quite complicated. 

Finally, another source of ambiguity concerns the choice of the basepoint xq. 
Winding the basepoint around some of the vortices (or some of the vortices around 
the basepoint) also changes our description of the multi-vortex configuration. This 
feature is not of much physical interest. When it is convenient, we will consider the 
basepoint to be far from all vortices, so that winding of vortices around the basepoint 
need not be considered. 

The gauge-invariant content of the classical n- vortex configuration can be encoded 
in the Wilson loop operators 

W^(^) (C) = -^x^''^ ^ exp ^ A • j j , (2.5) 

where x^'^^ denotes the character, and the dimension, of irreducible representation 
(z/). If C is a closed path around a single vortex, then W^'^\C), for all (z/), provides 
sufficient information to identify the class to which the vortex belongs. Likewise, 
W^W(C), for all (u) and all C in 7ri(Mfj) suffices to determine the n- vortex configura- 
tion, up to one overall conjugation. It is not sufficient though, to know W^'^\Ci) for 
all the generators Ci of 7ri(Mn); this determines only the class of each vortex, but not 
how the vortices are patched together. It is therefore more convenient and efficient to 
fix the gauge at the basepoint and assign group elements to standard paths around 
the vortices than to give a fully gauge-invariant description in terms of Wilson loops. 

When the unbroken gauge group H is discrete, one sometimes says that the theory 
respects a "local discrete symmetry." "'^'^ This terminology is used because a field 
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that is covariantly constant outside the vortex core is typically not single valued — on 
a closed path that encloses a vortex, a covariantly constant field is periodic only up 
to the action of an clement of H . (And, conversely, a field that is periodic cannot 
be covariantly constant.) But the phrase "local discrete symmetry" should not be 
misunderstood. Often, when we say that a symmetry is local, we mean that all 
physical states are required to be invariant under the symmetry. That is not what 
is meant here. Isolated H charges can exist, and can, in fact, be detected at infinite 
range via their Aharonov-Bohm interactions with vortices. 

2.2. Cheshire Charge 

Aside from their topological charges, vortices can also carry ordinary H quantum 
numbers!^'^^'"' These are easiest to discuss if we consider a vortex-antivortex pair, 
with trivial total topological charge. In this configuration, the vortex is described by 
a group element a, and the antivortex by a~^. But under H transformations, a mixes 
with other representatives of the conjugacy class to which it belongs. 

The action of H on the members of a class a defines a (reducible) representation 
which we denote as D^'^\ In D^'^\ each element of H is represented by a permutation 
of the class, according to 



This representation can be decomposed into irreducible representations of H. The 

physical interpretation is that the vortex pair can carry H charge. This charge is 
a property of the vortex-antivortex composite system and is not localized on either 
vortex or antivortex. It can be detected by means of the Aharonov-Bohm interaction 
of the composite with another vortex. 

For each class o; in if there is a unique state that can be constructed that is 
uncharged (transforms trivially under if); it is the superposition of group eigenstates 




a G a . 



(2.6) 




(2.7) 
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(where Ua denotes the order of the class) . Let us consider what happens to this state 
when an object that transforms as the irreducible representation (i/) of H passes 
between the two vortices. 

The composite system consisting of the vortex pair and the charged projectile has 
a well-defined if-charge; the composite transforms according to a representation of H 
that could be inferred by studying the Aharonov-Bohm interaction of the composite 
with other vortices. (It is to ensure this that we consider a pair of vortices with trivial 

total topological charge.) Interactions between the vortex pair and the projectile 
cannot change this total iJ-charge. From this observation, we can infer how the 
charge of the vortex pair changes when a charged object winds around one of the 

,. [30] 

vortices. 

It is convenient to imagine that the vortex pair actually interacts with a particle- 
antiparticle pair, combined in the trivial representation of H, so that the total H- 
charge of particles plus vortices is trivial. Then, if one particle winds around one vor- 
tex as in Fig. 2, the charge transferred to the vortex pair is opposite to the charge 
transferred to the particle pair; the vortices are in a state transforming as the repre- 
sentation {/J,)*, if the particle pair transforms as (/i). 

The charge-zero state of the particle-antiparticle pair has the group-theoretic 
structure 



(summed over i) where {e^ } denotes a basis for the vector space on which the 
representation (i/) acts, and rii, is the dimension of that representation. After the 
particle winds around a vortex with flux a & H, this state becomes 






(2.9) 



We may regard this state as a vector in the space on which the representation Z}(^)* (g) 
D^^^^ acts, expanded in terms of the natural basis ej'^-'* ® ej*^ /• This vector can also 
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be expanded in terms of the basis that block-diagonahzes the representation, as 

-^DM(a) = 5^c^,,^(^)'^. (2.10) 

Here, i; is a unit vector in the irreducible subspace acted on by the representation 
{n),y', the index y occurs because a given irreducible representation may occur in 
D^'^^ ® D^^^* more than once. 

We may interpret 

Pm = EI^/^>^|' (2-11) 

V 

as the probability that the final state of the particle pair transforms as D^^^^; corre- 
spondingly, it is the probability that, after interacting with the projectile, the vortex 
pair transforms as D^^^* (The |c|'s do not depend on how the class representative a 
is chosen.) Of course, can be nonvanishing only if is one of the irreducible 

constituents of D^^\ 

Thus we see that charge is transferred to the vortex pair, via the nonabelian 
Aharonov-Bohm effect, when a charged particle winds around the vortex; an 
initially uncharged vortex pair becomes a nontrivial superposition of charge eigen- 
states. This is "Cheshire charge." The charge is a global property of the vortex pair, 
and is not localized on the vortices. 



2.3. Strings 



The above description of vortices is easily translated into appropriate language 
to describe loops of string. A multi-string configuration can be characterized by 
assigning elements of H to each of a set of standard paths (beginning and ending at 
a basepoint xq) that link the string loops. Again, the paths that begin and end at 
xq and avoid the cores of all strings fall into homotopy classes; these classes are the 
elements of 7ri(M, xq), where M now denotes B? with all string loops removed. 
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Thus, when many strings arc present, there is an ambiguity in the choice of the 
"standard path" that hnks a given string, and a corresponding ambiguity in the 
assignment of group elements to strings. The physics associated with this ambiguity 
is that there is a long-range holonomy interaction between string loops. 

In Fig. 3a, standard paths a and (5 have been introduced that wind around two 
string loops. Suppose that, initially the homomorphism maps a to the group element 
a, and maps jS to the group element h. Now suppose that the h string loops winds 
through the a loop, in the sense defined by a. During the winding, the path shown in 
Fig. 3b is dragged to j3. Thus, the fiux associated with this path before the winding is 
the same as the flux associated with (3 after the winding. This path can be expressed 
as a(5a~^^ and so is mapped by the initial homomorphism to the group element aba~^ . 
We conclude that the final homomorphism after the winding maps (5 to 

h' = aha-^ , (2.12) 

and maps o; to a' = a. When a loop of h string winds through a loop of a string, its 
fiux becomes conjugated — it becomes an aba~^ string. 

For future reference, we remark that Fig. 3 has an alternative interpretation. The 
path shown in Fig. 3b will also be dragged to (3 if the a string winds around the 
basepoint xq, without the b loop passing through the a loop at all. Of course, the 
basepoint xq is completely arbitrary, and selected by convention, so this process has 
no actual physical effect on the b loop. Yet, the effect of parallel transport around 
the path P will be different than before, after the a loop has lassoed the basepoint. 
We will return to this point in Section 5. 

Furthermore, a string loop, like a pair of vortices, can carry an H charge, and can 
exchange charge with other charged objects by means of the Aharonov-Bohm effect. 
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2.4. Branching and entanglement 

Another property of strings will be relevant to our subsequent analysis; namely, 
that strings can branch. This feature is illustrated in Fig. 4. There a string labeled 
by c = a6 splits into two strings with labels a and b. If the gauge group is nonabelian, 
then the precise rules for composing the group elements at a branch point depend 
on our conventions — these (arbitrary) conventions associate each string with a path 
around the string that begins and ends at the basepoint, as shown in Fig. 4. 

Related to the branching phenomenon is another generic feature of nonabelian 
cosmic strings — strings that are labeled by noncommuting group elements become 

[31] 

entangled when they cross. Fig. 5 illustrates that when an a string and a b string 
cross, they become connected by a segment of c string. The mathematics underlying 
this entanglement is the very similar to that underlying the holonomy interaction 
between two nonabelian vortices, in two spatial dimensions. 

Let us calculate the flux c carried by the string segment that connects the a string 
and the b string, after the crossing. First, we establish our conventions by choosing 
standard paths a, /3, and 7 that encircle the string, as in Fig. 5a. The group element 
c associated with 7 can be determined if we demand that the bridge connecting the 
strings can be removed by re-crossing the strings. However, even once we have fixed 
all our conventions, the element c is not uniquely determined — there are two possible 
choices. This ambiguity arises because a crossing event involving two oriented strings 
has a handedness. A useful way to think about the handedness of the crossing is to 
imagine that the a and b strings are actually large closed loops. When the string 
loops cross, they become linked. But the linking number can be either +1 or — 1. 
(That is, we may consider a surface bounded by one of the loops. This surface inherits 
an orientation from the loop, defined (say) by the right-hand rule. If the other loop 
pierces the surface in the same sense as the outward-pointing normal, the linking 
number is otherwise, —1.) The two possible linking numbers correspond to the 
two possible types of crossing events. 

For one of the two types of crossings (the case in which the the loops have linking 
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number —1), the paths /? and /?' in Fig. 5b must be mapped to the same group 
element. In terms of the standard paths defined in Fig. 5a, we have [3' = a~^'jPa. 
(Recall our notation — in a composition of paths, the path on the right is traced first.) 
We therefore find b = a~^cba, or 

c = aba~^b~^ . (2.13) 

The connection with the holonomy interaction between vortices is clarified by Fig. 5c, 
where standard paths a^~^ and 7/3 are shown that wind around the "upper halves" 
of the entangled strings. Since ac~^ = {ab)a{ab)^^ and cb = {ab)b{ab)~'^ , the "flux" 
carried by the upper half strings differs from the flux carried by the lower half strings 
by conjugation by ab (just as winding an a vortex counterclockwise around a b vortex 
causes the flux of both vortices to be conjugated by ab). 

For the other type of crossing (the case in which the loops have linking number 
+1), the paths P and P" in Fig. 5d must be mapped to the same group element. This 
path is P" — aPja~^, so that b — abca~^, or 

c = b-^a~^ba . (2.14) 

Now the upper half strings carry flux (ab)~^a{ab) and ab~^b{ab), which is just the 
change in a two-vortex state that results from winding an a vortex around a b vortex 
in the clockwise sense. 
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3. Nonabelian walls bounded by strings 



If a gauge theory respects a local discrete symmetry G, then, as we have seen, 
the theory admits cosmic strings that arc classified by the elements of G. But if the 
group G is spontaneously broken to a subgroup H, then some of these strings become 
boundaries of domain walls. Here we will consider some of the properties of domain 
walls bounded by strings. It is necessary to understand these properties, in order 
to interpret the behavior of our correlation functions, and to use them to probe the 
phase structure of a model. 

In this section, as before, we assume that a standard choice of gauge has been 
made at a standard basepoint. This choice allows us to assign a definite group element 
to each cosmic string, and also fixes the embedding of the unbroken group if in G. 

3.1. Classifying walls bounded by string 

If a discrete gauge group G is spontaneously broken to a subgroup H, then the 
cosmic strings of the theory are classified by elements of H. The elements of G that 
are not in H are associated not with isolated strings, but rather with strings that 
bound segments of domain wall. 

There is not, however, a one-to-one correspondence between elements of G that 
are not in H and the physically distinguishable types of wall. This is because a single 
type of wall can be bounded by more than one type of string. As is illustrated in 
Fig. 6, strings labeled by a and b can bound the same wall, if h — ab~^ is in the 
unbroken group H. The composite of the two strings is an object that does not end 
on a wall, which means that the wall that ends on a a string must also end on a 6 
string (with conventions chosen as in the the figure). Thus, no locally measurable 
property distinguishes a wall that ends on an a string from a wall that ends on a 6 
string, if ab"^ G H. In other words, walls are classified by the cosets {gH} of in G 
(with the trivial coset corresponding to the trivial wall, i.e., no wall at all). Whether 
the walls are labeled by left or right cosets is a matter of convention. 
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Another way to think about this classification is in terms of the discontinuity of 

the order parameter across the domain wall. Suppose that the Higgs field $ that 
drives the breakdown from G to H transforms as the representation R of G, so 
that D^^^{h)^ = ^ ioi h E H. Across the wall labeled by a, the order parameter 
jumps from $ to D^^\a)^. This discontinuity is independent of the choice of coset 
representative. 

To avoid confusion, we should remark that the classification that we have just 
described is not the usual classification of domain walls in a spontaneously broken 
gauge theory. The point is that we are not classifying the domain walls that are 
absolutely stable, as in the standard analysis of topological defects in gauge theories. 
Instead, we are classifying the different walls that can end on strings. 

Notice that the criterion for two strings to be boundaries of the same wall is noi the 
same as the criterion for two strings to be (locally) indistinguishable objects. Strings 

labeled by a and h have the same structure if a and h are related by a = hbh~^ for 
some h & H. Two elements of G can be in the same H coset without being related in 
this way. Therefore, the same type of wall can have more than one type of boundary, 
in general. 

3.2. Wall decay 

The walls described above, those associated with elements of G that are not in 
the unbroken group H, are unstable. Such a wall will decay by nucleating a loop of 
string by quantum tunneling; the loop then expands, consuming the wall'^'. 

As we have seen, it is possible for one type of wall to end on more than one type of 
string. Hence, the wall may be able to decay in more than one way. Furthermore, if a 

large segment of wall bounded by an a string decays, two different decay modes that 
are locally indistinguishable may be globally distinguishable. For example, it may be 
that b and b' are both in the same H coset as a, and also that b' = hbh^^, for some 
h & H. Then the wall can decay by nucleating either a b string or a b' string, and 
the strings look the same locally. But ab~^ and ab'~^ need not be the same element 
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of H. Then the ribbon of wall that is produced by the nucleation of the loop (see 
Fig. 6), is different in the two cases. 

If a large loop of a string is prepared, which bounds a wall, then the wall can decay 
in any of the available channels. The string with the lowest tension is the most likely 
to nucleate, but all strings that can bound the wall nucleate at some nonvanishing 
rate. Many holes eventually appear in the wall, and the boundaries of the holes (the 
strings) expand due to the tension in the wall. Eventually the holes collide. 

When an a hole meets a b hole, junctions form, and an ab~^ string appears that 
bridges the hole. (See Fig. 7.) It may be that the b string has higher tension than the 
a string. Then the junctions will get pulled around the boundary of the hole. They 
annihilate, liberating a ab~^ string from the decaying wall system. 

This scenario shows that, when the wall bounding a very large loop of string 
decays, we may regard the decay process as dominated by the strings on which the 
wall can end that have the smallest possible tension. This insight will help us to 
interpret the behavior of the order parameter A^E, C). If an operator creates a world 
sheet of an a string, which bounds a wall, then the Aharonov-Bohm interaction with 
a charge whose world line winds around the string world sheet will be the same as the 
Aharonov-Bohm interaction of the charge with the ab~^ e H string, if the b string is 
the liglitest one that can nucleate on the wall. 

For completeness, we note another structure that can arise in walls bounded by 
strings. We have seen that an ab string can branch into an a string and a b string. It 
may be that all three strings are boundaries of walls. Then a possible configuration 
is shown in Fig. 8. Here the ab string is a vein in the wall, where the type of wall 
bounded by the a string joins the type of wall bounded by the b string. 
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4. The abelian order parameter 



The "abelian" order parameter introduced in Ref. 5 can be used to probe the 
reahzation of the center of the gauge group. Here we will briefly describe how this 
order parameter is used to distinguish the various phases of a Z^r gauge-spin system 
on the lattice. 

The model we consider has gauge variables Ui residing on the links (labeled by 
Z) of a cubic lattice, and spin variables (f)i residing on the sites (labeled hy i). Both 
gauge and spin variables take values in 

Zn = {exp{2TTik/N), k ^ 0,1,2, N - 1} , (4.1) 

The (Euclidean) action of the model is 

S — 'S'gauge + 'S'spin ) 

where 

5gauge = -/?J](?7p + C.C.) , (4.2) 
P 

and 

N-l 

^spin = - E E i[i^~'UcP)ir + C.C.) . (4.3) 

m=l I 

Here Up — YlieP^l associates with each elementary plaquette (labeled by P) the 
(ordered) product of the four Ufs associated with the (oriented) hnks of the plaquette, 
and ((f)~^U(f))ij — (f)~^Uij(f)j, for each pair ij of nearest-neighbor sites. In eq. (4.3), 
we have introduced an independent spin coupling constant corresponding to each 
nontrivial irreducible representation of Z^. 

Now we must consider how the operator -F(S) is to be defined in this lattice 
model.* Recall that inserting F(E) into a Green function is supposed to be equivalent 

* The analogous construction of the 't Hooft loop operator on the lattice was first discussed in 
Ref. 33 and Ref. 34. 
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to introducing a classical cosmic string source on the world sheet E. On the lattice 

(in 3+1 dimensions), we consider S to be a closed surface made up of plaquettes of 
the dual lattice. There is a set E* of plaquettes of the original lattice that are dual 
to the plaquettes of E.(See Fig. 9.) The operator Fn{'S) may be expressed as 

F^E) = n exp (/?(e2-WiVt;^ _ ^ ^ ^ ^^ (4 4) 

PeT.* 

In other words, to evaluate the path integral for a Green function with an insertion of 
F„(E), we modify the plaquette action on the plaquettes that are dual to E, according 
to 

Up e^^^'^/^C/p, P e E* . (4.5) 

This transformation of S'gaugc is equivalent to introducing n units of Zjsf magnetic 
flux on each of the plaquettes in E*. 

When the surface E is large, the vacuum expectation value of Fn behaves like 

(F^(E)) ~ exp (-«;^^")Area(E)) . (4.6) 

We may interpret k^'^™'* as the renormalization of the tension of the string source, due 
to quantum fluctuations of the charged matter flelds. Eq. (4.6) also has an alternative 
interpretation. We may think of E as a surface lying in a time slice, rather than as the 
world sheet of a string propagating through spacetime. Then Fn(E) is the operator 

F„(E) = exp (^^gs) . (4.7) 

where is the Z^v charge contained inside the surface E. Virtual pairs of charged 
particles near E cause the charge to fluctuate. If the theory has a mass gap, then 
the charge fluctuations near two elements of E become very weakly correlated when 
the elements are distantly separated. Thus, charge fluctuations cause the character- 
istic "area-law" falloff of (F„(E)) in eq. (4.6). 



21 



Consider the case P » 1 and « 1 (for all m). In this case, the gauge 
variables are highly ordered, and so we expect that charge is not confined. Fur- 
thermore, the spins are disordered, so there is no Higgs condensate to screen 
charge either. Thus, free Zjq charges should exist. We can check whether this expec- 
tation is correct by using perturbative expansions to evaluate (^(S, C)). 

Because the gauge coupling is weak, a frustrated plaquette (one with Up ^ 1) 
is very costly. An insertion of Fn(S) tends to frustrate the plaquettes in E*. These 
frustrations can be avoided, though, if the gauge variables assume a suitable con- 
figuration. Choose a three-dimensional hypersurface Q whose boundary is S. This 
hypersurface consists of a set of cubes of the dual lattice. Dual to these cubes is a 
set fl* of links of the original lattice. (See Fig. 9.) By choosing 

jj _ ^2nin/N ; g Q* 

(4.8) 

Ui= 1 , l^Q\ 

we can avoid frustrating any plaquettes. (This is the lattice analog of a "singular 
gauge transformation" that removes the string "singularity" on E.) With the hnks 
in the configuration eq. (4.8) the fundamental representation Wilson loop operator 

W(C) ^l[Ui (4.9) 
leC 

acquires the value exp{2'Kink / N) , where k is the linking number of C and E.(See 
Fig. 10.) 

Of course, this choice of link variables changes the nearest neighbor interactions 
for each pair of spins that is joined by a link contained in fl*. But the spins are strongly 
coupled and highly disordered, so that they are nearly indifferent to this change; it 
is much more costly to frustrate a plaquette than to frustrate a link. Therefore, the 
expectation value of the ABOP 

is dominated by small fluctuations about the "background" eq. (4.8), and we conclude 
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that (^A^^\j:,C)j satisfies cq. (1.2). Thus, there is an infinite range Aharonov-Bohm 
interaction, and free charges exist. (Further details of this analysis may be found 
in Ref. 5.) 

If the gauge couphng is strong (/3 << 1), then there are no free Zj^ charges. The 
Wilson loop introduces a Z^ charge as a classical source, but confinement causes a 
pair of Z]\f charges to be produced, so that the charge is dynamicaUy shielded. We 
expect that eq. (1.3) is satisfied, and this can be verified in the small-/? expansion'^'. 

It is also interesting to consider the case in which P is large and some of the 
spin couplings ^ym are also large. Then the matter field is ordered and the matter 
"condensate" screens the charge. 

For example, suppose that 

7m' >> 1 , m' = m , 

(4.11) 

7m' << 1 , m . 

In effect, then, the operator (f/^ condenses, and Zjsi is broken down to the kernel of 
the representation (m); this is Zm-, where M is the greatest common factor of N 
and m. We anticipate, therefore, that the operator -Fn(S) introduces a stable cosmic 
string world sheet provided that nm = [modN) (so that the flux carried by the 
string is in the unbroken group Zm)- Otherwise, the string introduced by Fji(E) is 
the boundary of a domain wall. This wall is unstable, and decays by nucleating a 
loop of string. 

//we assume that the flux of the string created by Fn{Ti) combines with the flux 
of the nucleated string to give a trivial total flux, then we anticipate that the ABOP 
behaves as 

lim {aI:\e, C)) = exp (^-^k{^, C)^ , mn = (mod N) , 
lim (A„(E, C)) = 1 , mn ^ (mod N) . 

From this behavior, we could easily infer that the unbroken symmetry is Zm] a string 
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has a nontrivial Aharonov-Bohm interaction with a charge if and only if the flux 
carried by the string is in Zm- 

However, as we emphasized in Section 3.3, even if the string introduced by F„ 
is the boundary of a domain wall, and the wall decays by nucleating a string loop, 
it need not be the case that the nucleated string combined with the classical string 
source has trivial flux. Thus, eq. (4.12) does not, in general, correctly describe the 
behavior of the ABOP in the hmit eq. (4.11). To see what actually happens, let us 
analyze the consequences of eq. (4.11) using perturbative expansions. 

For 7^ >> 1, spins with nontrivial Zm charge are highly ordered, and frustrating 
these spins is very costly. Now there is a competition between the reluctance of the 
system to frustrate a plaquette (when (3 » 1), and its reluctance to frustrate a Zm 
spin. If the operator F„(E) is inserted, then, as we have seen, frustrated plaquettes 
can be avoided if the links of Q* are excited. But, if exp{27rin/N) ^ Zm-, we can excite 
the Ui^ on these links only at the cost of frustrating the spins there. The number of 
links in Q* increases like the volume enclosed by E, and the number of plaquettes in 
E* increases only as the area of E. Thus, for E sufficiently large, frustrated plaquettes 
are favored over frustrated links. This is just a realization on the lattice of the decay 
of a domain wall by nucleation of a loop of string, where the frustrated spins comprise 
the wall, and the frustrated plaquettes comprise the nucleated string. 

But as we discussed in Section 3.3, it is sometimes possible for a domain wall 
to decay in more than one way. To illustrate the possibilities, we will consider two 

different special cases. 

4.1. Z^^Z2 

Consider, for example, a Z4 model. For /9 >> 1 and 72 >> 1, the Z4 symmetry 
will break to Z2. One finds, indeed, that l^J^^^ iJ^^C^ behaves as in eq. (4.12). 

Understanding the behavior of J^^^ and J^^^ involves a subtlety. The n = 3 
string is the anti-string of the n = 1 string; therefore, they both have the same 
tension. Furthermore, exp(27ri/4) and exp(— 27ri/4) belong to the same coset of Z2 
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in Z4, so the n = 1 and n = 3 strings are both boundaries of the same domain wall. 
When the operator -Fi(S) is inserted, the resulting domain wall can decay in two 
different ways. One way, the composite string that is created has trivial Z2 flux; the 
other way, it has nontrivial Z2 flux. 

Correspondingly, when -Fi(S) is inserted, we may choose either Ui — 1 or Ui — —1 
on the links of O*; the weak spin coupling in the action (the 72 term) depends only on 
U^, and so is not frustrated either way. But wc also need to consider the dependence 
on the strong couplings 71 and 73. Expanding in these small parameters, one finds 
that the effective tension of the composite string is rcnormalized by spin fiuctuations. 
The renormalization raises the tension of the composite string with nontrivial Z2 flux 
relative to that with trivial flux. Thus, the conflguration such that the composite 
string has trivial Z2 flux really does dominate when -Fi(E) (or Fs{T,)) is inserted. 
Therefore, A^^ and A^^ do behave as in eq. (4.12). 

4.2. Zq^Zs 

Now consider a Zq model. For P » 1 and 73 >> 1, the Zq symmetry will break 
to Z3. ^42^^^ and A^^'^ behave as in eq. (4.12); so do A^^^ and a''^\ 

Now consider a'^^ . The n — 3 string that is introduced by F3 is the boundary of 
a domain wall. But this wall can terminate on an n = 1 or n = 5 string, as well as 
on an n = 3 string. Furthermore, the n — 1 (or n — 5) string has lower tension than 
the n — 3 string, so that nucleation of this string is favored. 

In other words, the conflgurations that dominate, when ^3(2) is inserted, have 
Ui — exp(27ri/3) (or Ui — exp(— 27ri/3)) for Z e Q* By exciting the links on Q*, we 
reduce the degree of frustration of the plaquettes of E. And we do so at negligible 
cost, because the m — 3 term in 5'spin cannot distinguish Ui — exp(±27ri/3) from 

Thus, the nucleated string, combined with the classical string source, has non- 
trivial flux. The combined flux can be either n — 2 or n — A, and both occur with 
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equal probability. And therefore, even though F3 inserts a string that decays, A3 
shows nontrivial behavior; we have 



lim (^a'^s\^^ = l (f''''''^ + e-^'^*'^/^^ - cos(27ri//3) 



(4.13) 



(if E and C have linking number A; = 1). 

Thus, for dynamical reasons, eq. (4.12) is not satisfied. Nevertheless, the behavior 

of An^ has an unambiguous interpretation. For example, the v- dependence of (a^o^) 



shows that the effective string introduced by F3 has flux that takes values in the Z3 
subgroup of Zq. The only possible interpretation is that the unbroken subgroup is 



Indeed, the behavior of i^An j always contains sufficient information to unam- 
biguously identify the unbroken subgroup of an abelian discrete local symmetry group. 

Other order parameters have been suggested that probe the realization of the cen- 

[IT 20 18] 

ter of the gauge group. ' ' We will comment on the efficacy of these in Section 9.3. 



We have seen in the previous section how the Aharonov-Bohm effect can be 
used to probe the phase structure of an abelian gauge theory. We now want to 
extend this procedure, so that it can be used in nonabelian theories. The basic 
strategy will be the same as before — we will introduce strings and charges as classical 
sources, and investigate the dynamical response of the theory to these sources. But 
the implementation of this strategy is more delicate in the nonabelian case. The main 
stumbling block is the problem of introducing nonabelian cosmic strings in a lattice 
gauge theory, which we address in this section. 

We will also construct operators that create or annihilate dynamical string loops; 
correlation functions of these can be used to study the dynamical properties of strings. 




^3. 




5. Nonabelian strings on the lattice 
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5.1. String calibration 

We consider a theory with (discrete) gauge group G. We are interested in how 
the local G symmetry is realized. Specifically, we wish to identify the subgroup H 
of G that admits free charges. (G quantum numbers may be confined, or may be 
screened by a Higgs condensate.) 

We investigate the realization of G by assembling a laboratory that is equipped 
with cosmic string loops. As described in Ref. 21 we can cahbrate the string loops 
with a beam of particles that transform as some faithful (not necessarily irreducible) 
representation {K) of G. We choose an arbitrary basepoint xq, and a basis for the 
representation (i?) at that point. We direct the beam from the basepoint to a beam 
splitter, allow the two beams to pass on either side of the string, and then recombine 
the beams and study the resultant interference pattern. 

If the string is in a "group eigenstate" with fiux (as defined in eq. (2.1)) a E G, and 
\u) is the wave- function in internal-symmetry space of a particle at the basepoint, 
then, when the particle is transported around a closed path that begins and ends at 
xq, the wave-function is modified according to 

\u)^D^^\a) \u) . (5.1) 

By observing the interference pattern, we can measure 

{u\D^^\a)\u) . (5.2) 

By varying \u), we can then determine all matrix elements of D^^\a), and hence a 
itself. Given our choice of basepoint, and a choice of a basis for D^^\ this procedure 
allows us to associate a well-defined group element with each string. 

In the case where free G charges may not exist, we must be careful that the sep- 
aration between the two beams remains small compared to any confinement distance 
scale or Higgs screening length. At the same time, of course, the separation must be 
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large compared to the thickness of the string core; we assume that the core is small 
compared to the characteristic distance scale of the dynamics that we wish to study. 

Once we have calibrated the strings by measuring their Aharonov-Bohm inter- 
actions with nearby charges, we probe the dynamics of the theory by measuring the 
Aharonov-Bohm interactions of strings with distant charges. In this way, we hope to 
learn what quantum numbers are confined or screened, and to infer the "unbroken" 
subgroup H. 

(As we have noted, each string loop associated with an element of G that is not in 
H will become the boundary of a domain wall. Thus it might seem that a good way 
to identify H is to observe which strings are attached to walls. Indeed, at sufficiently 
weak coupling, this is a sensible procedure, because the walls are very long lived. But 
at intermediate coupling this procedure may fail, because the walls decay rapidly by 
nucleating string loops. The time scale for the decay may be comparable to the time 
required to assemble and calibrate a string loop.) 

Some problems with the calibration procedure should be pointed out. The first 
is that a string loop associated with a definite group element is not an eigenstate of 
the Hamiltonian of the theory. An a string and a b string will mix with each other 

[21] 

if a and b are in the same H conjugacy class, and the energy eigenstates will be 
linear combinations of group eigenstates that transform as definite representations of 
the unbroken group H. (We will have more to say about this mixing in Section 7.2.) 
But the time scale for this mixing increases exponentially with the size of the string 
loop. For our purposes, it will usually be legitimate to ignore the mixing and regard 
the strings as group eigenstates. 

But there is another more serious problem. While the quantum fluctuations that 

change the identity of a string loop are very rare (when the loop is large), there are 
other, much less rare, quantum fluctuations that can change the Aharonov-Bohm 
phase that is detected. Suppose, for example, that a virtual b loop nucleates, lassoes 
the basepoint xq and then reannihilates. Naturally, this process has no physical 
effect on an a loop that is far from the basepoint. Yet it changes the Aharonov- 
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Bohm phase acquired by a particle in representation (z/) that winds around the string, 
beginning and ending at xq; wc saw in Section 2.3 that the phase becomes D^^\hab~^) 
rather than Z^*^^) (a). (This process is depicted, in 2+1 dimensions, in Fig. 11.) Such 
fluctuations are suppressed at sufficiently weak coupling, but they are present, at 
some level, for any finite coupling. They result in an ambiguity in the Aharonov- 
Bohm phase associated with a string, even if we fix the basepoint and a basis for the 
faithful representation {v). 

To avoid this ambiguity, we are forced to take a trace of the representation; the 
character x^^\'A is unaffected by these fluctuations. Thus, the existence of free G 
charges (and of an unscreened Aharonov-Bohm interaction) can be probed by the 
Wilson loop operator eq. (2.5). (And since W^^\C) is gauge invariant, there is no 
need for the loop C to contain the basepoint xq.) Let Fa(E,xo) denote the operator 
that inserts a string world sheet on E; the string is associated with a & G, relative 
to the basepoint xq. (We consider below how this operator is constructed.) Then we 
may define the ABOP 

Aa (S'-o,C)-^^^^^^^^^^^^(^^^^^^ . (5.3) 

In a phase with free G charges, we expect that 

hm (4"^(E,a;o;C)) = — x(^)(a'=(^'^)) , (5.4) 



n 



V 



where k denotes the linking number. (Even this statement requires a qualification, 
for the loop C must not be permitted to come close to retracing itself on successive 
passages around E.) 

We will consider in Section 9 how ^^'^^(E, xq; C) is expected to behave for other 
realizations of the local G symmetry. 
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5.2. Inserting string world sheets 

Now we will consider in more detail the problem of constructing an operator that 
inserts nonabelian strings on specified world sheets. The construction will be guided 
by the discussion in Section 2 of the general formalism, and by the discussion above 
of the calibration procedure. 

Given a set of disjoint closed two-dimensional surfaces Ei, E2 . . ., this operator 
is to introduce string world sheets on these surfaces, where the strings are associated 
with the group elements 01,02,.... As we have emphasized, the strings must be 
referred to a common basepoint xq, and each associated group element depends on 
the choice of a standard path that begins and ends at xq and winds around the 
string. The group elements oi, 02, . . . are defined up to one overall conjugation, which 
corresponds to a gauge transformation at the basepoint xq. 

Furthermore, the definition of our operator must be insensitive to the possible 
confinement or screening of G quantum numbers; it must specify the short-distance 
structure of the string, and leave it up to the dynamics of the theory whether the 
string can be detected at long range. It must realize the calibration procedure de- 
scribed above, in which the short-range Aharonov-Bohm interactions of the string 
are determined. 

We wish to define an operator Faj^a2,...,a„(^i) ^2, ■ ■ ■ , ^ru 3;o) so that an inser- 
tion of -Fai,02,...,0n(Si) ^2, ■ ■ ■ , xq) in a Green function introduces strings on world 
sheets Ei, E2, . . . , E^j. We will first describe the operator in heuristic terms, and 
then give a more precise description in the context of lattice gauge theory. Loosely 
speaking, an insertion of -Fai,a2,...,a„(Si) ^2, ■ ■ ■ , xq) in a Green function imposes 
a restriction on the gauge field configurations that are included in the path integral. 
Make a choice of a standard path that begins at xq, winds around the surface Ei, 
and returns to xq. Now consider a path Pi homotopic to the standard path that runs 
from Xq to a point on Ei, traverses an infinitesimal closed loop that encloses Ei, and 
then retraces itself, returning to xq. When Fai,02,...,o„(Si) ^2, ■ ■ ■ , ^n, xq) is inserted. 
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the gauge field configuration is restricted to satisfy 



P exp 



i j 



\ 



A ■ dx 



(5.5) 



for any such Pi. Similar restrictions apply for paths P2, Ps, . . . , Pn that wind around 
E2, E3, . . . , E„, as in Fig. 12. These restrictions on the path integral define the oper- 
ator Fa^^a2,-.-,a„ 

The operator so constructed is not gauge invariant. Of course, if it is inserted 
in a Green function with gauge-invariant operators, its gauge-invariant part will be 
projected out. Alternatively we may obtain an explicitly gauge-invariant operator by 
averaging over gauge transformations at the basepoint xq, obtaining 



— ^ Fgaig-\ga2g-\...,gang-^ (Ei, E2, . . . , E^, Xq) , 



(5.6) 



where no is the order of the group. 
We may now consider the operator 



Aai,a2,...,a„0-'l, '^2, • • • , Xq] C) 



i^ai,a,,...,a„(Si,E2,...,E„Xo)iy('')(C) 



(Ei,E2,...,E„,xo))(iyM(C)) 



(5.7) 



In a phase with free G charge, if the loop C is homotopic to PiPjP^ . . ., we have 



lim (Ai^Ja2,...,a„(Si,E2, . . . ,E„,xo;C)) = —x^^^iaiajak - . .) . 

Til/ 



(5.8) 



Thus, when a charge passes around several string loops in succession, the Aharonov- 
Bohm phases acquired in each successive passage are combined coherently. The co- 
herence is maintained because we have defined the various loops in reference to the 
same basepoint xq. If different basepoints had been chosen instead, then the group 
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elements associated with the various string loops would have been averaged over 
conjugacy classes independently. Since 



-^^DM(^a^-i) = -^X^^)(a)l (5.9) 

g€G 

(which follows from Schur's lemma), we find, for example, 

(Foi(Ei,xo)Fa2(E2,yo))(W^^'^-'(C^)) 

if the loop C winds around world sheets Ei and E2 in succession. 

1) ^o) ^T^d (E2,|/o) have distinct basepoints {xq ^ yo), then each by 
itself introduces a gaugc-singlet object. Inserting both operators combines two string 
loops as (trivial) charge eigenstatcs, rather than as group eigenstates. This is the 
reason for the lack of coherence in the Aharonov-Bohm interaction characterized by 
eq. (5.10). To combine two strings that are not both gauge singlets, we must include 
in our operator some nonlocal construct that bridges the gap between the string loops, 
just as a string of electric flux must be included in a gauge-invariant operator that 
creates a widely separated quark-antiquark pair. In the case of two (or more) string 
loops, this nonlocal connection between the loops is provided by referring the loops 
to a common basepoint. 

It is important, actually, that the Aharonov-Bohm interactions of an additional 
loop combine incoherently with the Aharonov-Bohm interactions of existing loops, 
when the group element associated with the additional loop is averaged over a con- 
jugacy class. It is this incoherence property that ensures that the effects of virtual 
string loops do not spoil the "factorization up to a phase" represented by eq. (5.8). 

We will now describe the construction of -Fa(S,a;o) in more detail, by specifying 
how the construction is carried out in a lattice theory with discrete gauge group G. As 
discussed in Section 4, E is to be regarded as a closed surface consisting of plaquettes 
of the dual lattice; these plaquettes are dual to a set E* of plaquettes of the original 
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lattice. For each plaquctte P e E*, we choose a path Ip on the lattice that connects 
the bascpoint xq to one of the corners of P. These paths are chosen so that each 
closed loop IpPlp^ is "homotopic" to the standard loop that links S. The various 
paths may be chosen arbitrarily, except that the union of all the paths should not 
contain any closed loops. 

The effect of the operator Fa(T,,xo) is to modify the plaquette action on each 
plaquette in E*. Suppose, for example, that the plaquette action is 

sSLp--Px^''\Up) + c.c. (5.11) 

(where i? is a representation of G that must be specified to define the theory). Then 
an insertion of Fa{T^,xo) modifies the action according to 

4auge,P -^^^""^ {yi.<^Vi;'Up) + C.C., P e E*, (5.12) 

where 

Vi, = llUi. (5.13) 

lelp 

Alternatively, we may write 

Fa(E, xq) = n (/^^^"^^ (Vlp(^Vi;'Up^ - Px^'^HUp) + C.C.) . (5.14) 
PeE* 

The operator Faj^a2,...,o„(Si) ^2, ■ ■ ■ , E^^, xq) that inserts many string loops is con- 
structed by a straightforward generalization of this procedure. 

As constructed, Pa(E, xq) is not gauge invariant. When inserted in a Green func- 
tion with gauge-invariant operators, though, it has the same effect as the explicitly 
gauge-invariant operator in which a is averaged over a conjugacy class, as in eq. (5.6). 
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It is also instructive to consider the correlator of Fa(E,a;o) with the operator 



where the product is taken over a closed set of links that begins at ends at xq. The 
trace of U'^''\C,xo) is {riu times) the Wilson loop operator W^^^C). But U^^''\C,xo) 
itself, like Fa{12, xq), is not invariant under a gauge transformation that acts nontriv- 
ially at the basepoint xq. 

In a phase with free G charges, and in the leading order of weak coupling pertur- 
bation theory, one finds that* 



This equation merely states that, once a loop of string has been calibrated, the same 
Aharonov-Bohm phase can be recovered again if another interference experiment is 
subsequently performed. But we have already emphasized that quantum fluctuations 
(such as the virtual string depicted in Fig. 11) can spoil this result. Indeed, when 
higher orders in the weak coupling expansion are included, it is seen that the correlator 
of Fo(E,xo) and U^'^\C,xq) fails to "factorize" as in eq. (5.16), even when E and C 
are far apart. As stated before, we must consider the correlator of Fa(S, xq) with the 
gauge-invariant operator W^'^\C), in order to extract an Aharonov-Bohm "phase" 
that depends only on the topological linking of E and C (in the limit of infinite 
separation) . 

To conclude this section, we must ask how the operator Fa(E, xq) depends on the 
basepoint xq and on the choice of the paths {Ip}- (We consider correlation functions 
of Fa with gauge-invariant operators, so we regard Fa as gauge-invariant, with a 
averaged over a conjugacy class.) We note that, since the union Up Ip of all the 

★ This calculation is described in more detail in Section 6.2 and in the Appendix. 




(5.15) 




(5.16) 
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paths contains no closed loops, we can choose a gauge with Ui = 1 for all I e Ip. In 
this gauge, the group clement a is inserted directly on the plaquettes of S*. Thus, 
it is clear that the correlator of Fa(E,a;o) with any local gauge-invariant operator is 
independent of the choice of basepoint and paths. 

The correlator of Fai,02,...,o„(Si) ^2, ■ ■ ■ , ^ru ^^o) with nonlocal gauge-invariant op- 
erators (hke Wilson loop operators) depends on the choice of path only to the extent 
that we have already noted. That is, it depends on the choice of the "standard paths" 
that enter the calibration of the string loops.We may change the paths from the base- 
point xq to the plaquettes of Eg, by winding these paths around Si, as in Fig. 13. 
In effect, this change alters the group element assigned to the string world sheet on 
E2; 02 becomes replaced by aia2a^^. As we noted in Section 2, this ambiguity in 
(El, E2, . . . , Xq) is the origin of the holonomy interaction between string 
loops (or vortices) 

5.3. Inserting string loops 

An insertion of the operator Fa(E,a;o) introduces a string on the closed world 
sheet E. The string may be regarded as an infinitely heavy classical source. Thus 
F(j(E,xo) is closely analogous to the Wilson loop operator W^'^\C)] an insertion of 
W^'^\C) introduces an infinitely heavy classical source (in representation (v)) on the 
closed world line C. 

But we will also find use at times for an operator that creates (or annihilates) a 
dynamical cosmic string. Such an operator can be obtained by a simple modification of 

the construction described above — namely, the surface E is chosen to be, rather than 
a closed surface, a surface with nontrivial boundary C.(See Fig. 14.) Carrying out 
the same procedure as before for such a surface, we arrive at an operator Ba{C, E, xq) 
that creates a string (or annihilates an anti-string) on the loop C. More generally, an 
operator 

-Bai,a2,...,an(Cl, ^i, C2, E2, . . . , C^, E^, Xq) (5-17) 

creates strings on the loops Ci, C2, . . . , C„, with all strings referred to a common 
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basepoint xq. 

This construction generalizes a construction devised by 't Hooftj^*^' and we will 
refer to B as the '"t Hooft loop" operator. However, 't Hooft considered the situation 
in which the strings have no Aharonov-Bohm interactions with other fields. In that 
case, the surface E is an invisible gauge artifact, and B{C) depends on C alone. We 
are interested in strings that have Aharonov-Bohm interactions, and in that case E 
is not invisible.'^' 

To better understand why the 't Hooft operator Ba{C, E, xq) must depend on the 
surface E as well as on the loop C, it is helpful to think about a theory defined in 2+1 
spacetime dimensions. In that case, the operator Ba{x, y, E, xq) creates a vortex at x 
and an anti- vortex at y; E is a path connecting x and y. But no gauge-invariant local 
operator exists that creates an isolated vortex at x. Because a vortex can be detected 
at infinite range via the Aharonov-Bohm effect, there is a vortex superselection rule. 
Hence, the operator that creates a vortex cannot be local; it has a semi-infinite string 
that can be seen by the fields of the theory. (Similarly, there is an electric charge 
superselection rule in quantum electrodynamics. No gauge-invariant local operator 
can create an isolated electron; the operator that creates an electron must also create 
a string of electric fiux that ends on the electron.) 

Correlation functions such as (i?aj^a2(Ci, Ei, C2, S2, a;o)) can be used to deter- 
mine the tension of a dynamical string, or the amplitude for mixing between group 
eigenstates, as we will describe in Section 7. 

(As an aside, we remark that an 't Hooft loop, like a Wilson loop, admits an 
alternative interpretation. If a Wilson loop operator acts on a timelike slice, it is 
natural to interpret it as an insertion of a classical charged source, as noted above. 
But if the Wilson loop acts on a spacclike slice, we may interpret it as an operator 
that creates a closed electric flux tube. We have noted that we may think of an 
't Hooft operator acting on a spacelike slice as an object that creates a cosmic string. 
Alternatively, we may interpret the an 't Hooft loop acting on a timelike slice as an 
insertion of a classical magnetic monopole source. In the situation originally consid- 



36 



ered by 't Hooft, the monopole satisfied the Dirac quantization condition, and so its 
"Dirac string" was invisible. We are considering a situation in which the Dirac string 
is visible; the surface E bounded by C is the world sheet of this Dirac string.) 



6. Classical strings in the pure gauge theory 



We will now analyze the behavior of the operator A^a (S, xq; C), using perturba- 
tive expansions. Here we will consider the case of a pure gauge theory with (discrete) 
gauge group G. In Section 9, we will consider the effects of introducing matter. 

6.1. Strong coupling 

Although our main interest is in the physics at weak gauge coupling, we will make 
a few comments about the strong-coupling behavior of the pure gauge theory. 

The plaquette action of the theory is taken to be 



where is a representation of G. For /3 << 1, this theory confines sources that 
transform as certain irreducible representations of G. The criterion for a source to be 
confined is easiest to state in the case where R is irreducible. In that case, a source in 
the irreducible representation (z/) is not confined if and only if there are non-negative 

integers ki and k2 such that 



where R* denotes the complex conjugate of i?, and (1) denotes the trivial represen- 
tation. The point is that if cq. (6.2) is satisfied, then it is possible for a source in 
representation {y) to be "screened by ghions." In other words, the electric flux tube 
that terminates on a {v) source can break due to glue fluctuations. (The criterion 



gauge,P 



(6.1) 



(i?)^i ® (i?*)^2 D (1) , 



(6.2) 
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eq. (6.2) generalizes the familiar notion that an adjoint representation source is un- 

confined in a strongly coupled SU{N) gauge theory.) If (z/) is confined, then the 
expectation value of the corresponding Wilson loop operator decays for a large loop 
C like 

(iy(^)(C)) ~exp(-«;MArea(C)) , (6.3) 

where Area(C) is the minimal area of a surface bounded by C, and is the tension 
of a flux tube that carries electric flux in the representation (z/). (Of course, virtual 
ghic may partially screen the source; k^'^^ is the tension of the lightest flux tube that 
can terminate on a (z/) source.) 

Since the gauge variables can see a classical string with "magnetic flux" a e G, 
the quantum fluctuations of the gauge variables renormalize the tension of the string 
by an amount Ka^^"^ ■ Thus, the expectation value of the operator F decays for a large 
surface E hke 

(F,(E, xo)) ~ exp (-4''"^ A(E)) , (6.4) 

where A(E) is the area of E. The calculation of the leading contribution to K,a^^\ for 
/3 << 1, is described in the Appendix. 

Now consider the operator A^\Tj,xq;C). Since there arc no infinite-range 
Aharonov-Bohm interactions in the confining phase of the theory, we might expect 
that 

lim(Ai")(S,a;o;C)) = l. (6.5) 

This is indeed found for a representation (u) that is not confined (such that the electric 
fiux tube can break). Different behavior is found, however, for a representation (z/) 
that is confined. The operator A^a\^, xq; C) can have a nontrivial expectation value 
because the electric fiux tube stretched across C crosses the surface E at certain 
points. In the leading order of strong-coupling perturbation theory (and assuming 
that no "partial screening" of the source occurs), each crossing contributes to A the 
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factor 




(6.6) 



(or its complex conjugate, depending on the relative orientation of the flux tube and 
S at the point of crossing). Thus, even when S and C are far apart, {a[I^\j1, xq; C)) 
is not a purely topological quantity that depends only on the linking number of S 
and C. Of course, in higher orders in the strong coupling expansion, the behavior of 
A^a^ becomes still more complicated. 

Similarly, if a loop C links with two different surfaces Si and E2, the operator 

A^ai,a2^l^ S2,a;o; C) acquires a factor {1 / ny)x^'^\o'i) each time the electric flux tube 
crosses Si and a factor / ni/)x^^\ci2) each time it crosses S2 (in the leading order 
of the strong-coupling expansion). Because of confinement, the string world sheets 
combine as trivial charge eigenstates rather than group eigenstates, even though both 
are defined with respect to a common basepoint. 

6.2. Weak coupling 

For 13 » 1 there is no confinement, and the Wilson loop operator exhibits 
perimeter law decay for any representation {u). 

The operator Fa(E,xo) decays as in eq. (6.4). The calculation of the leading 
behavior of K^a^^\ for /3 >> 1, is described in the Appendix. 

When the operator Fa(S,a;o) is inserted, there is a configuration of the gauge 
variables such that no plaquettes are excited. This configuration can be constructed 
by choosing a set VL of cubes of the dual lattice such that the boundary of VL is E. 
Each cube in Q is dual to a hnk of the original lattice. The configuration with no 
excited plaquettes is 



This configuration is unique up to a gauge transformation. (The gauge transforma- 
tions are deformations of O.) 



Ui^a , I eQ* , 



(6.7) 
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Weak-coupling perturbation theory is carried out by expanding in the number 
of excited plaquettes, and in the degree of excitation. In the hmit /3 — > cxo, the con- 
figurations with the minimal number of excited plaquettes dominate. By calculating 
the Wilson loop operator for the configuration with no excited plaquettes, we verify 
eq. (5.4) in the weak-coupling limit. Thus we find, as expected, that G charges are 
neither confined nor screened. Similarly, we may verify eq. (5.8) and eq. (5.16) in this 
limit. 

When higher-order corrections in weak-coupling perturbation theory are com- 
puted, we find as anticipated that eq. (5.4) and eq. (5.8) continue to hold. But 
eq. (5.16) does not survive. These corrections are further discussed in the Appendix. 

We wish to make one other remark here about the weak-coupling expansion, 
which might help to avoid confusion. To calculate the weak-coupling behav- 
ior of correlation functions that involve the operator -Faj^aj, .•■•an(^i? ^2, • • • , Sn,a;o), 
we first construct the configuration that has no excited plaquettes when 
-Poi,a2,...,an(^l) ^2, • • • , xq) is inserted. This construction is a straightforward gen- 
eralization of eq. (6. 7). However, there are topologically inequivalent ways of choosing 
nonintersecting surfaces f^i, . . . , that are bounded by Si, . . . , E^j (as in Fig 15). 
Thus, one might get the impression that there can be two (or more) gauge-inequivalent 
configurations that both have no excited plaquettes. But this is not the case. To see 
why not, it is important to keep track of the basepoint, and of the paths from the 
basepoint to the loops. If Qi is distorted past Q2, as in Fig. 15, then Qi crosses the 
paths from the basepoint to E2. To avoid exciting any plaquettes, then, the links 
contained in must now take the value 01020^^. (This, again, is a reflection of the 
holonomy interaction between string loops.) Therefore, a Wilson loop that crosses 
fli first and Q2 second, in Fig. 15a, behaves exactly the same way as a Wilson loop 
that crosses Q2 first and Qi second, in Fig. 15b. There is a unique gauge equivalence 
class of configurations with no excited plaquettes, just as there should be. 
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7. Dynamical strings and vortices 



We will now consider how the 't Hooft operator can be used to investigate the 
properties of dynamical stnngs (in 3+1 dimensions) and vortices (in 2+1 dimensions). 

7.1. String tension and vortex mass 

The operator 5^ ,^-1 (Ci, Ei, C2, S2, xq) can be used to compute the tension of 
a cosmic string that carries magnetic flux a (with a caveat described below). This 
operator creates an a string on Ci and annihilates it on C2. Thus, when the loops 
are large and far apart, we have 

(^a,a-i(Ci,Si,C2,S2,a;o)) 

- exp (^-^i^r\A{Ei) + ^(E2))) exp (-/.^^'^^^(Ci, C2)) . ^^'^^ 

Here is the renormalization of the tension of a "classical" string source, and 

4^^^^) is the tension of a dynamical string; A{Ci, C2) is the area of the minimal surface 
with boundary Ci U C2.(See Fig. 16.) If the loops Ci and C2 are chosen to be far 
apart compared to the correlation length of the theory, but close together compared 
to the size of the loops, then the dependence of {Ba^a~^{Ci,Ti\,C2-,^2-,XQ)) on the 
separation between the loops determines the tension n'a^'^^ ■ 

Actually, the same information can be extracted from the behavior of the simpler 
operator Ba{C, E, xq). For a large loop C, we have 

(S«(C, E, xo)) ~ exp {-Kt'"'^ A{i:)) exp (-nf^'^^ A{C)) , (7.2) 

where A{C) is the area of the minimal surface bounded by C.(See Fig. 17.) Since A(E) 
and A{C) can be varied independently, tih^"^^ can be determined. (Or, Fa(S,a;o) can 
be used to measure K"a^'^\) The calculation of in weak-coupling perturbation 

theory is described in the Appendix. (In the strong-coupling limit, we have k^'^^^^ — 
— there are no stable magnetic flux tubes.) 
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Obviously, the same procedure can be used to calculate the mass of a vortex, in 
2+1 dimensions. 

The existence of a stable string (or vortex) can itself be used to probe the phase 
structure of the theory. In a confining phase, stable magnetic flux tubes do not exist; 
they "melt" due to magnetic disorder. If the G gauge symmetry is spontaneously 
broken to a subgroup i7, stable a strings exist only if a e if. Otherwise, an a string 
is the boundary of a domain wall, which decays as described in Section 3.3. 

If no stable a string exists, then Ba{C, E, xq) might not create a stable string. If 
it does not, its expectation value will behave, for a large loop C, like 

{Ba{C, S, xo)) ~ exp exp (-m^'^)p(C)) , (7.3) 

where P{C) is the perimeter of C. It may seem, then, that by measuring 
{Ba{C,Ti,XQ)), and determining whether it decays as in eq. (7.2) or as in eq. (7.3), 
we can find out whether a is contained in the unbroken subgroup H or not. However, 
there are subtleties. One problem is that there may be tradeoff between the depen- 
dence of (Sa(C, E, xo)) on A(E) and its dependence on A{C). Eq. (7.3) will apply if 
the domain wall bounded by E decays by nucleating an string that completely 
cancels the flux of the classical string source on E. But it may be that K^a^^^ can be 
reduced if the nucleated string only partially screens the flux of the source. (We saw 
an instance of this phenomenon in the Zq example that was discussed in Section 4.) 
The advantage gained from reducing /t^'^™) may more than compensate for the cost 
of a nonvanishing n^'^^'^^] then Ba will decay as in eq. (7.2), even though a ^ H. 

Another complication can arise if the unbroken subgroup H is not a normal 
subgroup of G. For then a typical G conjugacy class contains both elements that are 
in H and elements that are not in H. Recall that i?a(C, E,xo) is actually averaged 
over the G conjugacy class that contains a. One particular iJ-class contained in 
this G-class will dominate the asymptotic behavior of {BaiCj'EjXo)), and whether 
eq. (7.3) or eq. (7.2) applies depends on which class dominates. 
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We will return to the problem of finding a suitable order parameter, that can be 
used to identify H, in Section 9. 



7.2. Mixing 

If the discrete gauge group G is unbroken, then elements of G in the same con- 
jugacy class are associated with strings that transform into each other under the 
action of G. In the classical limit, these "group eigenstate" strings are degenerate 
energy eigcnstates. (There is also a further degeneracy associated with "parity," 
which changes the orientation, and so transforms the a string into the string.) 
Quantum mechanically, these states mix with one another, and the degeneracy is 
lifted. The true energy eigenstates are "charge eigenstates" that transform according 

[21] 

to irreducible representations of G (and parity). 

This mixing can be computed using the 't Hooft operator. We consider a corre- 
lation function in which an a loop is created on Ci and a b^^ loop is annihilated on 
C2, where b = gag~^ for some g E G. It is crucial that the two strings be defined 
with respect to the same basepoint xq. Otherwise, we would average a and b over the 
conjugacy class independently, and the correlation function would be dominated by 
the propagation of an a string from Ci to C2, rather than the mixing of an a string 
with a b string. 

Let Ci and C2 be two congruent loops, one directly above the other as inFig. 18. 
The separation between the loops is large compared to the correlation length of the 
theory, but small compared to the size of the loops. If the ab~^ string is stable, 
then, in the weak coupling limit, the correlation function will be dominated by the 
configuration in Fig 18a. In this configuration, the world sheets of the a and b strings 
join, and the loop at which they join is the boundary of the world sheet of an ab~^ 
string. If this configuration dominates, then 



Ba,b{Gi, El, C2, E2, xo) ~ exp (-«f"")^(Ei) - 4''"^A(E2)) 



(7.4) 
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But if the ab~^ string is unstable, and decays to a widely separated a string and 
string, then the configuration in Fig. 18b will dominate. Here, the world sheets of 
the dynamical strings are stretched tightly across Ci and C2. If this configuration 
dominates, then 



exp(-(«:?^'^) + «:f'^V(Ci)) 



(7.5) 



Thus, the mixing amplitude in the weak-coupling limit is either 

exp (-(4'^") + 4'^^'^V(^^i)) (7.6) 



o(mix) 

g a,b rsj 



or 



(^(mix) 



e^P ^(Ci)) , (7.7) 



whichever is larger. 

Actually, these are not the most general possibilities, for the ab~^ string may be 
unstable, and may prefer to decay in some other channel. If the ab^^ string decays 
to n widely separated stable strings that carry flux ci, C2, . . . , c„, then we flnd 



e-«ir^exp(-(^4*"')A(C,)) . (7.8) 

These results were previously derived in Ref. 21. 
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8. Interactions 



8.1. Entanglement and holonomy interactions 

The Aharonov-Bohm interaction between a classical string source and a classical 
charged source was studied in Section 6. Correlation functions can also be used to 
study the holonomy interaction between two vortices (in 2+1 dimensions), or two 
strings (in 3+1 dimensions), and also the entanglement of strings. 

In 2+1 dimensions, two loops Ei and E2 can link, as shown inFig. 19. Then 
the quantity (F^ ^(Si, S2, xq)) probes what happens when one vortex source winds 
around another. As we noted in Section 2, there is a holonomy interaction between 
the vortices if a and b do not commute; both are conjugated by afc'^"''^*'. This means 
that the vortex world lines cannot close (the vortex pairs cannot re-annihilate) unless 
there is an exchange of topological quantum numbers between the two vortices. 

An {ab)a{ab)^^ vortex becomes an a vortex, and an {ab)b{ab)^^ vortex becomes 
a b vortex, if the flux aba~^b~^ flows from the b vortex to the a vortex* Therefore, 
the quantum numbers that must be exchanged are those of an aba~^b~^ vortex. We 
thus find, in the weak coupling limit, 

(-^a,6(Sl,S2,Xo) 

~ exp (-mi^^")p(Ei) - m^r^PiJ:^)) exp {-m^X-M^i, E2)) ■ ^^'^^ 

Here 

^(ren) 

the renormalization of the mass of a classical vortex source, m^'^^^^ is 
the mass of a dynamical vortex, P(E) is the perimeter of E, and I/(Ei, E2) is the 
length of the shortest path that connects Si and S2. (If the aba^^b^^ vortex is 
unstable, then ^^fco-ife-i replaced by the sum of the masses of the vortices to which 
it decays.) 

★ Here we have used different conventions to assign group elements to the entangled vortex world 
lines than the conventions used in Section 2.4 to assign group elements to entangled strings in 
three spatial dimensions. 
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The leading behavior of weak-couphng perturbation theory on the lattice is found 
by identifying the configurations with the minimal number of excited plaquettes, as 
we described in Section 6. To find the leading contribution to {Fa^b{^l:^2, xq)) , we 
choose surfaces fli and fl2 that are bounded by Ei and S2 respectively, and construct 
a configuration with Ui — a on links in 0| and Ui — b on links in Og- But where fli 
and intersect, there is a string of excited plaquettes with 

42ge,p--/5x^''^(«^'«"'^'"') + c-c. (8.2) 

(See Fig. 20.) By choosing fli and fl2 so that this string has minimal length, we 
obtain eq. (8.1). 

(If the aba~^b~^ vortex is unstable, then the dominant configuration in the weak- 
coupling limit is found by a slightly modified procedure. Either Ei or E2, or both, 
becomes the boundary of several surfaces or 0,2- The fink configuration is chosen 
so that Ui — Ci on Q^* and Ui — dj on Qg'^, where Yi^i — a and Yi^j — b. These 
surfaces intersect along several strings that connect the two world lines; thus, the 
classical vortices exchange several separated dynamical vortices.) 

In 3+1 dimensions, strings labeled by noncommuting group elements become 
entangled when they cross, as we described in Section 2.4. String world sheets in 
four dimensions generically intersect at isolated points. (A point of intersection is 
a type of "instanton." ) Consider surfaces Ei and E2 that cross at two points, as 
shown inFig. 21. Because the strings entangle, an a world sheet on Ei and a b world 
sheet on E2 become joined by the world sheet of an aba~^b~^ string. Thus, in the 
weak-coupling limit, we find 

{Fa,bi^l,^2,Xo) 

~ exp [-K^r^AiE,) - K^r^AiE^)) exp (-4^-6-^(^1' ^2)) , ^^'^^ 

where A(Ei, E2) is the area of the minimal surface that joins Ei and E2. Again, this 
result is easily verified using weak-coupling perturbation theory on the lattice. 
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It is also instructive to consider (Fa^;,(Ei, E2, xq)) where the surfaces Ei and E2 
have the topology shown in Fig. 22. Here Ei is a torus that links once with the 
sphere E2. As before, Qi and Q2 unavoidably intersect along a surface of frustrated 
plaquettes; we find the leading behavior of (^^ ^(Ei, E2,a:o)) by choosing fli and fl2 
so that this intersection has minimal area. The result is again eq. (8.3). But now, if 
the width of the torus is small compared to the size of the sphere, ^(Ei, E2) is the 
area of the minimal slice through the torus. 

To interpret this result, we recall the observation in Section 2.3, that when a loop 
of a string winds around a loop of b string, it becomes a loop of bab^^ string. We may 
regard Fig. 22 as a depiction of a process in spacetime, in which loops of a and 
string are produced, the a loop winds around the b string, and the a and strings 
then annihilate. But because of the Aharonov-Bohm interaction between the strings, 
this process is disallowed; an loop cannot annihilate a bab~^ loop. Therefore, 
(Fa^6(Ei, E2, xq)) is suppressed by the amplitude for the bab~^ string to become, via 
quantum tunneling, an a string. Indeed, comparing with eq. (7.7), we see that 

e-^r;:- ^ exp (-4^-6-^^""^ (Si)) (8-4) 
is the suppression factor in eq. (8.3). 
8.2. Cheshire charge and the Borromean rings 

As discussed in Section 2.2, it is an inevitable consequence of the nonabelian 
Aharonov-Bohm phenomenon that a loop of string (or a pair of vortices) can carry 

[5 12 13] 

charge, and can exchange charge with other charged objects. ' ' . Let us see how 
this property is reflected in the behavior of the correlation functions of our string 
operators.'*"' 

Consider the following process (in 2+1 dimensions), depicted in Fig. 23. First, 
a pair consisting of an a vortex and its anti-vortex spontaneously appears; the total 
charge of the pair is trivial. Then a charged particle in the representation {u) winds 
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counterclockwise around the a vortex; thus, charge is transferred to the vortex pair, 
as described in Section 2.2. Next, a b vortex winds around the (charged) vortex pair, 
acquiring an Aharonov-Bohm phase that is sensitive to the charge of the pair. Then 
the {v) particle winds clockwise around the a vortex, discharging the pair. Finally, 
the a vortex and anti- vortex reunite and annihilate. 

If the vortices and charged particle are treated as classical sources, this process 
is captured by the correlation function (^F^ ;,(Ei, E2, xq) W^'^\C)'^, where the three 
loops El, E2 and C are configured as in Fig. 23. This is a topologically nontrivial 
joining of three loops known as the "Borromean rings;" no two loops are linked, yet 
the loops cannot be separated without crossing. Because this correlation function 
has the spacetime interpretation described above, we may anticipate that, when the 
three loops are joined but far apart, ;,(Ei, E2, xq) W^'^\C)^ differs from its value 
for large unjoined loops by a topological multiplicative factor — the Aharonov-Bohm 
phase acquired by the b vortex that winds around the charged vortex pair. 

To calculate ^^'^^(Ei, E2, xq; C), we proceed, as usual, by finding the configuration 
of the link variables such that, when Fa^b is inserted, there are no frustrated plaquettes. 
We pick surfaces fli, bounded by Ei, and fl2, bounded by E2, that do not intersect; 
such surfaces are shown in Fig. 24. We then choose (up to a gauge transformation) 

Ui^a , / e fii , 

Ui = b, leQ*2, (8.5) 

Now we compute W^'^\C) in this configuration. As is clear in Fig 24, the loop C 
crosses first Qi in a positive sense, then Q2 in a negative sense, then Qi in a negative 
sense, and finally Q2 in a positive sense, before closing. The corresponding path- 
ordered exponential is ba~^b~^a, and taking a trace yields'^"' 

hm (a[^1{Ei,E2,xo;C)) ~ — x^'Haba-^b-') . (8.6) 

\ ' / Til/ 

We should verify that the factor eq. (8.6) can be interpreted as the Aharonov- 
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Bohm phase acquired by a 6 vortex that winds around a vortex pair with Cheshire 
charge. The interpretation is easiest if we exphcitly average a over the representatives 
of the class to which it belongs. If this averaging is not performed, then there is a 
correlation between the choice of the class representative that is used to measure the 
charge, and the choice of the class representative for the string that is being measured; 
this correlation makes the interpretation of the measurement more complicated. Of 
course, averaging a and h over class representatives independently is equivalent to 
defining a and h with reference to distinct basepoints; we have 



(Fa(Ei,a;o)) (Ffe(S2,2/o)) {W^^liC)) no 

(8.7) 

We recall from the discussion in Section 2.2 that, when a particle in representation 
(i/) winds around one of the vortices of a pair that is initially uncharged, the final 
state of the vortex pair is a superposition of states with various values of the charge. 
Eq. (8.7) gives the expectation value of the Aharonov-Bohm phase acquired by a 6 
vortex that winds around the charged pair. If is the probability that the vortex 
pair transforms as the irreducible representation /x, then we have 

E^^M -X^'\b) = — E - X^''\gag-'bga-'g-'b-') . (8.8) 

Indeed, eq.(8.8) can be verified directly, withp^ given by eq. (2.11) and (2.10). Details 
will be presented elsewhere.'^"' 

The above discussion can be applied, with hardly any modification, to the case 
of charged string loops, in 3+1 dimensions. There is an analog of the Borromean 
ring configuration, in which two disjoint closed surfaces are joined by a closed loop, 
although the loop is not hnked with either surface. Eq. (8.7) holds, if Ei^2 are surfaces, 
and C is a loop, in this configuration. 
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9. Gauge-Higgs systems 



9.1. Weak coupling 

In this section, we will consider models with Higgs fields coupled to gauge fields. 
We wish to investigate the screening of gauge charges — and of Aharonov-Bohm 
interactions — due to the Higgs mechanism. 

For a gauge system with (discrete) gauge group G, we introduce a Higgs field (pi, 
defined on sites, that takes values in G. We then include in the action of the model 
the term 

^Higgs = -E>E (^^"^^ +C.c) , (9.1) 

M I 

where the sum runs over all irreducible representations (//) of G. (Compare cq. (4.3).) 
Suppose that some of the 7^'s arc large (>> 1), while all others arc small (or zero). 
Then, we can analyze the behavior of this model using perturbative methods. 

In the weak-coupling limit 7^ — > 00, 0~^C/0 becomes restricted to the kernel of 
the representation (//) at each hnk, and so C/p is also "frozen" to the kernel at each 
plaquette. Thus, G becomes "spontaneously broken" to the subgroup i/ = Ker(i^M). 
By choosing several of the 7^'s to be large, we can break G to the kernel of a reducible 
representation. Indeed, breakdown to any normal subgroup of G can be achieved in 
this way, since every normal subgroup is the kernel of some representation. 

Of course, we know from the usual continuum weak-coupling analysis that more 
general patterns of symmetry breakdown are possible. In that analysis, the Higgs 
field resides in the vector space on which a representation of G acts, and, in principle, 
the unbroken subgroup could be the stability group of any vector in this space. In the 
model with S'jjiggs given by eq. (9.1), these more general patterns are not obtained 
when all of the Higgs couplings cither vanish or arc very large. They might, of course 
be obtained at intermediate coupling. After suitable block-spin transformations are 
performed, the effective theory that describes the infrared behavior of the model 
would be a "continuum" theory for which the usual analysis could apply. 
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Anyway, with the breaking of G to a normal subgroup H implemented as de- 
scribed above, we can proceed to calculate the ABOP A^a'^iT,, xq; C) in weak-coupling 
perturbation theory, and so probe the fate of the Aharonov-Bohm interaction in the 
Higgs model. The calculation gives the expected results. But there is one difficulty. 

The problem is that the irreducible representation (z/) of G is typically a re- 
ducible representation of the subgroup H. In general, just one of the irreducible 
representations of H that is contained in (z/) will dominate the asymptotic behavior 
of 1^(^)(C) when the loop C is large. So it is the Aharonov-Bohm interaction of this 
H representation with an a string that is probed by A^a"^ . In fact, if {v) contains the 
trivial H representation, then this will dominate at weak coupling, and ^^1^*^^^ will 
behave trivially even though {v) may contain other H representations that do have 
Aharonov-Bohm interactions with the string. 

We have encountered a general problem with the interpretation of I^A^a ''^ that 
arises whenever G is nonabelian and is partially broken. If I^A^a^'^ behaves nontrivially 
for some choice of a, then we know that an unscreened charge must dominate W^^^] 
by varying a, we can obtain information about the representation according to which 
this charge transforms. But if l^A^a^'^ behaves trivially for all a, we know only that 
W^^^ is dominated by a screened charge. We have no a priori knowledge of how this 
charge transforms, or of how other H representations contained in (z/) interact with 
strings. This problem complicates the task of inferring H from the behavior of our 
correlation functions. 

9.2. String stability 

We noted in Section 7.1 that, if G is a discrete group, then the existence (or 
not) of a stable string with flux a E G provides a criterion for determining whether 
a is contained in the unbroken group H. It may seem, then, that the operator 
Ba{G, E, xq), which creates an a string, can be used to probe the realization of the G 
gauge symmetry. But there are two problems. 
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First, when inserted in a Green function with gauge-invariant operators, Ba is 
actually averaged over the G conjugacy class to which a belongs. This class may 
contain some elements that are in the unbroken group H and some that are not, 
which complicates the interpretation of {Ba). 

Second, even if no element of the class that contains a is in H, Ba{C, E, xq) may 
nevertheless create a stable string. This can happen, as we saw in Section 4.2, due 
to the competition between the renormalization of the tension of the classical string 

that propagates on S and the tension of the dynamical string whose world sheet is 
bounded by C. 

Let us ignore the second problem for the moment, and address the first. We first 
note that the situation is relatively simple if if is a normal subgroup of G. In that 
case, a given G conjugacy class is either entirely contained in H or is disjoint from 
H. By studying Ba for various class representatives, we can determine which classes 
are contained in H, and so reconstruct H. 

The general case is somewhat more complicated, and Ba by itself is not sufficient 
to completely determine the unbroken subgroup H. Instead, one way to proceed is 
the following: If G = {ai, 02, ... , Onc} ^ finite group of order no, consider 



in which all group elements appear. Choose the loops to be large and widely separated, 
with their sizes ordered according to 



Now Soi,02,... is effectively averaged over gauge transformations at the basepoint, and 
so may be replaced by 



{Bai,a2,...,a„Q{Cl, Ei, C2, E2, . . . , C, 




(9.2) 




(9.3) 




'gaig ^,5025 ^,...,ganQg 



i(Ci,Ei,C2,Il2,...,CnG,IlfiQ,Xo) . 



(9.4) 



52 



The expectation value is dominated by the configurations that minimize 

+ (9.8) 

(where k^^^^ — for a ^ H). This condition may not determine g uniquely. 
Among those g that minimize eq. (9.5), the dominant configurations are such that 
K,a2^^A{T,2) + K,^^'^^ A{C2) is also minimized. And so on. Now by varying A{C) and 
A{T,) independently, we find the group elements a for which k-^^J^i — 0. We thus 
determine the unbroken subgroup H up to one overall conjugacy H — > gHg~^. This 
ambiguity is expected; it corresponds to the freedom to change the embedding of H 
in G by performing a gauge transformation. (Of course, the calculation of eq. (9.2) 
involves an average with respect to this embedding.) 

In passing, we have found the tension of all of the stable strings associated with 
the various elements of H. 

If G is abelian, as in the discussion in Section 4.2, then we can overcome the 
second problem (that Ba may create a stable string even for a ^ G) easily enough. 
By measuring A^"^ , we can determine the Aharonov-Bohm interactions of the string 
created by Ba, and so identify the fiux of the string as belonging to H. If G is 
nonabelian, though, life is more complicated. For as we noted above, we have no a 
priori knowledge of what H representation dominates the asymptotic decay of W^'^\ 
While it seems altogether physically reasonable that the stabihty and Aharonov- 
Bohm interactions of strings can be used to identify an unbroken gauge group H, it 
is not so easily to specify how this should be done with gauge-invariant correlation 
functions. 
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9.3. The VOOP 



Another promising probe of charge screening in a gauge theory is the "vacuum 

[17 20 18] 

overlap order parameter" (VOOP) proposed by Marcu and Fredenhagen. ' ' Let 
us compare and contrast the VOOP with the Aharonov-Bohm order parameter that 
has been discussed in this paper. 

Suppose that an H gauge theory contains a matter field ^^'^^ that transforms as 
the irreducible representation (/i) of H. If the local H symmetry is unbroken, and 

the representation [ji) is not confined or screened, then (loosely speaking) the field 
<l>^ should create a stable particle. We can express this in gauge-invariant language. 
If X and y are distantly separated points, and Px,y is a path from y to x, consider the 
nonlocal gauge-invariant operator 



K^^^\x,y,Px,y) = 4^)t^(/^) n ■ (9-6) 



If $(A')t creates a free charge, then this charge must propagate between x and y, as 
inFig. 25. Thus, we have 

(x('^)(x,y,P))~exp(-Ma L(P)) exp (-M(;;i \x - y\) , (9.7) 

where M^^^ is the mass of the stable particle created by <|)(^)t^ M^e^ is the 
renormalization of the mass of the classical source propagating along P. (Here, L[P) 
is the length of P, and \x — y\ is the distance from x to y.) Since Mrel can be 
determined independently by measuring (^W^^\C)) (or by varying L(P) with — y| 



fixed), eq. (9.7) can be used to find M^y^- But if the representation (//) is screened (or 
confined), then (^K^^\x, y, P)^ will become independent of \x—y\ for large separation; 
in effect, we have M^^]^ = 0, with M^^^ defined by eq. (9.7). Thus, Marcu and 
Fredenhagen suggested that M^^^ > if and only if $(/^)t creates a free charge. 
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This construction is obviously closely similar to our discussion in Section 7.1 of 
string stability, and it suffers from related problems. Suppose that, in a model with 
gauge group G, spontaneous breakdown to a subgroup H occurs. Then, an irreducible 
representation (/x) of G contains various irreducible representations of H. It may 
be that among these representations are some that can exist as free charges, and 
others that are screened. We might expect, then, that the screened charges dominate 
is:(^)^,sothatM^5;^ = even though (//) contains some unscreened representations 
of H. Even this is not clear; because of the competition between Mdyn and Mren in 
eq. (9.7), the free charges may actually dominate. 

Because of these difficulties, it is not at all easy, in general, to identify the unbro- 
ken gauge group H based on the behavior of (^K^^\x, y, P)"^- 

9.4. Order parameters: Some concluding remarks 

The existence of stable cosmic strings (or vortices), and of Aharonov-Bohm in- 
teractions between strings and free charges, can be used to identify and classify the 
various phases of a gauge theory. Yet, because of the problems discussed above, it 
proves difficult to formulate a general procedure that unambiguously specifies the 
realization of the gauge symmetry, e.g., the "unbroken" subgroup. This is surprising 
(to us), but we are reluctant to attach any fundamental significance to it. 

In fact, though, just because of these problems, the phase structure of certain 
gauge theories may be richer than one might naively expect. For example, if the 
gauge group G is "spontaneously broken" to H, then, as we have remarked, just 
one of the irreducible H representations contained in the G representation (u) will 
dominate the asymptotic behavior of (w^'^\C)^. On some surface in the parameter 
space of the theory, a "crossover" may occur, where this H representation changes. 
Thus, the order parameter A^^ might be nonanalytic on this surface — the surface 
would be a sort of phase boundary, even though no change in symmetry occurs there. 

For most of the discussion in this paper, we have taken the gauge group G to be a 
discrete group. There is no obstacle, however, to generahzing our results to the case 
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where G is a continuous group. 



10. Electric flux tubes and dynamical monopoles 

In this paper, we have systematically discussed the Aharonov-Bohm interac- 
tions between magnetic flux tubes and electric charges, which can occur in the Higgs 
phase of a gauge theory. Central to the discussion has been the effect of quantum- 
mechanical electric charge fluctuations on the interaction. 

There is another type of Aharonov-Bohm interaction, which can occur in the 
confining phase of a gauge theory — the interaction between an electric flux tube and 
a magnetic charge. The existence of such interactions has been noted previouslyj^^'^^'^'^' 

[5 9] 

as have the implications concerning the "magnetic hair" carried by black holes. ' But 
the effects of quantum-mechanical magnetic charge fluctuations on this incarnation 
of the Aharonov-Bohm phenomenon have not been analyzed before. 

In this section, we will develop a quantum field-theoretic treatment of electric 
flux tubes in a confining gauge theory that contains dynamical magnetic monopoles, 
and will investigate the interactions of flux tubes with monopoles. This treatment, of 
course, will be closely similar to our theory of magnetic flux tubes, in a Higgs phase. 

First, we will describe (following closely the work of Srednicki and Susskind ) 
how dynamical magnetic monopoles are introduced into a lattice gauge theory. Then 
we will construct a generalized Wilson loop operator that creates an electric flux 
tube that has Aharonov-Bohm interactions with the monopoles. By studying the 
properties of the correlation functions of this operator, we will investigate the effect 
of magnetic charge fluctuations on the Aharonov-Bohm effect. 
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10.1. Magnetic Monopoles on the Lattice 

For definiteness, we will consider a gauge theory with gauge group SU{N). A 
pure SU{N) gauge theory (no matter), or a theory with the matter fields transform- 
ing trivially under the center Zjv of SU{N), admits magnetic monopoles with Zjv 
magnetic charges. (Fields that transform under Zn would be able to see the Dirac 
string of such a monopole; that is, the monopole with minimal Zj^ charge would 
not satisfy the Dirac quantization condition.) If SU{N) /Zjy is actually the unbro- 
ken gauge symmetry of an underlying theory with simply connected gauge group 
G D SU{N)/Zn, where G is broken via the Higgs mechanism, then such monopoles 
arise as topological solitons. For example, there are models with G — Spin{N'^ — 1) 
broken to SU{N)/Zjq that contain Zjv monopoles. 

The (usual) 't Hooft loop operator B{C) inserts a world line of a Zjv monopole 
along the closed loop C. But the monopole introduced by an 't Hooft loop is a classical 
source, not a dynamical object. We wish to elevate the monopoles to the status of 
dynamical degrees of freedom, and introduce a coupling constant that controls the 
effects of virtual monopoles. 

It is most convenient to choose the dynamical variable to be a "Dirac string field" 
that is summed over in the path integral. On the lattice, this field associates with 
each plaquette P a quantity rjp G Zjv, which may be regarded as the Zjsi magnetic 
flux carried by a Dirac string that pierces that plaquette. If the total magnetic flux 
entering a cube of the lattice is nontrivial, then a magnetic monopole resides in that 
cube. Of course, since the Dirac strings themselves must be unobservable, this theory 
should respect, as well as the usual SU{N) local symmetry, an additional Z^v local 
symmetry that deforms the Dirac strings (without, of course, changing where the 
monopoles are). 

[231 

The construction of this theory was described by Srednicki and Susskind. In 
the absence of matter flelds, it has the action 

S^-(3j2^r{r)pUp)-XY, Vc + cc. , (10.1) 

plaq cubes 
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where 

Pec 

is the product of the six Typ's associated with the oriented faces of the cube c; in other 
words, Tjc is the Zf^ magnetic charge inside c. 

The extra Zn local symmetry respected by this theory is defined on links; it acts 
according to 

(10.3) 

P3l, 

where n = 0,1,2, N — 1. This transformation moves the Dirac strings without 
changing the magnetic charge r]c or the magnetic flux rjpUp that appears in the 
gauge field plaquette action. (The quantity Up is not invariant under the extra local 
symmetry, and so is unphysical; it can be interpreted as a fictitious magnetic fiux 
that includes a contribution from the (unobservable) Dirac string that crosses P.) 

The coupling constant A controls the strength of the effects of virtual monopoles. 
For A << 1, magnetic charge fluctuations occur copiously, but for A >> 1, magnetic 
charge fluctuations are strongly suppressed. In the limit A — >■ oo, the monopoles 
freeze out, and eq. (10.1) becomes the usual Wilson action. 

Matter fields can be coupled to the gauge theory in the usual way. For example 
if (f)i e SU (N) is defined on sites (labeled by i) , we may define 

^matter = "7 X^""^ {{<P~'U<P)i) + C.C. . (10.4) 
links 

This is invariant under the local symmetry eq. (10.3) only if Z^ is contained in 
the kernel of the representation (R). Of course, this is just the Dirac quantization 
condition — the matter fields must be chosen so that the string of a monopole is 
invisible. We may introduce matter fields that are only invariant under some subgroup 
of Zn, but then we must restrict the rjp^s to take values in that subgroup. 
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10.2. Wilson Loop Operator 



Now we want to define a Wilson loop operator W^'^\C) that creates an electric 
flux tube on the loop C. But the usual construction 




(10.5) 



is not invariant under the local symmetry eq. (10.3) unless is contained in the ker- 
nel of the representation (u). The operator W^'^\C) is physically sensible, then, only 
if the string created by it has no Aharonov-Bohm interaction with a Zp^ monopole. 
In fact, in the SU{N) theory defined by eq. (10.1), this operator does not create a 
stable fiux tube at all — even at strong coupling (/3 << 1), glue fluctuations cause the 
string to break. 

An operator that creates a stable electric flux tube that does have Aharonov- 
Bohm interactions with monopoles cannot depend on the loop C alone; it must also 
depend on a surface S that is bounded by C. After our discussion of the 't Hooft 
operator in Section 5.3, this comes as no surprise. A (naive) Wilson loop operator, 
in the presence of monopoles, is a multi-valued object, for it acquires a non-trivial 
Aharonov-Bohm phase when the loop winds around a magnetic charge. To construct 
a single- valued object, we introduce a branch cut on the surface E, so that the operator 
jumps discontinuously when a monopole crosses E. We augment the naive Wilson 
loop operator, then, by a factor that counts the total magnetic flux of the Dirac 
strings that cross E, obtaining 

W^'^\C,E)^-X^-')((jlUi)[llvp)) . (10.6) 

This operator is invariant under eq. (10.3), for any representation {u) of SU{N), and 
of course it reduces to W^'^\C) if {v) represents Zjv trivially. 
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We may also consider the degenerate case in which the loop C shrinks to a point. 
If (i/) represents faithfully, then the operator 




(10.7) 



inserts the world sheet of a tube with minimal electric flux, as a classical source, 
on the surface S. (Of course, if (z/) represents Z^ trivially, then the flux tube is 
invisible, and the operator is trivial.) If we think of the surface E as lying in a time 
slice, then G(X1) has another interpretation; it is a magnetic charge operator that 
measures the total Z^ magnetic flux through E. Obviously, G(E) is the magnetic 
analog of the operator that we called -F(E) in Section 4. 

Note also that if matter is introduced as in cq. (10.4), an operator can be con- 
structed that creates a separated particle-antiparticle pair; it is 



where Pij is a path from the site i to the site j. This operator is invariant under the 
local symmetry, if (R) represents Z^ trivially. 

10.3. 't Hooft Loop Operator 

We can also construct an 't Hooft loop operator Bn{C); it inserts on the world 
line C a classical monopole source with Zjv magnetic charge n — 1,2, . . . , N — 1. 
Alternatively, we may interpret Bn{C), acting in a time slice, as an operator that 
creates a loop of magnetic flux tube on C. 

The construction of this 't Hooft operator may be carried out in much the same 
way as in a gauge theory without dynamical monopoles. Wc regard C as a closed loop 
composed of links of the dual lattice, and we chose an arbitrary surface E (composed 
of plaquettes of the dual lattice) whose boundary is C. Dual to the plaquettes of E is a 




(10.8) 
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set E* of plaquettes of the original lattice. Now, if there are no dynamical monopoles, 
we regard the ?7p's as classical (non-fiuctuating) variables. Then, to evaluate a Green 
function with an insertion of Bn{C), we modify the plaquette action on the plaquettes 
in E*, by choosing 

(10.9) 

When there are dynamical monopoles, however, and the r/p's are in the configuration 
eq. (10.9), the cubes that are dual to the links of C are frustrated. Thus, this 
configuration represents a (i?/riam«ca/ monopolc propagating on the world line C. Since 
we want the operator Bn{C) to introduce a classical monopole source, we should 
modify the cube action so that the configuration eq. (10.9) docs not frustrate any 
cubes. Thus, we evaluate a Green function with an insertion of Bn{C) by changing 
the cube action according to 

rj^ ^ e^^i^/^Tjc , c e C* . (10.10) 

Equivalently, we have 

BniC) = n '^^P (A(e2'^'"/^rye - 77c + cc)) . (10.11) 

cec* 

This operator is the magnetic analog of the Wilson loop operator. 

In the weakly-coupled gauge theory without dynamical monopoles, the 't Hooft 
loop operator creates a stable magnetic flux tube, and {BniC)) exhibits area-law 
decay. But when there are dynamical monopoles, Bn{C) always exhibits perimeter- 
law decay. The interpretation is clear. For any finite A, a Zjs[ magnetic flux tube is 
unstable, for the tube can break via nucleation of a monopole-antimonopole pair. 

Note that the operator defined by eq. (10.11) makes sense even if C is an open 
path rather than a closed loop. That is, we may define in like fashion an operator 

Bn{Pi,j)= n exp(A(e2-WiV,7e-77c + c.c)) , (10.12) 
where i and j are sites of the dual lattice, and Pij is a path connecting these sites. 
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This operator creates a monopole-antimonopolc pair, connected by a Dirac string; it 
can be used to compute the mass of a dynamical monopole. Obviously, it is closely 
analogous to the operator eq. (10.8). 

Having now in hand the operator Bn{C) that introduces a classical monopole 
on the world line C, and the operator G^'^\T) that introduces a classical Zjq electric 
flux tube on the world sheet E, we are ready to construct the operator, analogous 
to An\'^.,C), that probes the Aharonov-Bohm interaction between monopoles and 
electric flux tubes; it is 

pWh^ n - "^'"'(S) Bn(c) 

<='^)-(GM(E)> (B„(Cr 

If there is an infinite range Aharonov-Bohm interaction, the expectation value of this 
operator will have the asymptotic behavior 



hm {E'n'} = —x^"^ [ (e^«"/^^j ) , (10.14) 



where A;(E, C) is the linking number of E and C. 

Finally, we remark that it is straightforward to generalize the construction in 

Section 5, and define an operator that creates a loop of nonabelian cosmic string, in a 
theory that contains dynamical monopoles. We need only be cognizant of the change 
in the plaquette action that occurs when dynamical monopoles are included; a string 
that carries flux a e SU {N) is now created by 

S«(C,E,xo) = n exp|/3(tr(V^^ay^; Vt/p) -tr(77pC/p)) +C.C.} , (10.15) 

where E is a surface on the dual lattice, bounded by the loop C. When C shrinks to 
a point, we obtain the operator Fa(E,a;o) that introduces a classical string source on 
the closed world sheet E. 
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10.4. MoNOPOLE Condensation 

A pure SU (N) gauge theory, without dynamical monopoles, is confining at strong 
coupUng. For /3 << 1, gauge field fiuctuations are unsuppressed, and the resulting 
magnetic disorder gives rise to stable electric flux tubes. We want to explore how 
dynamical magnetic monopoles modify the physics of this theory. 

First, we consider the parameter regime P « 1 and A << 1, so that virtual 
monopoles are unsuppressed. It is easy to anticipate what will happen. A "monopole 
condensate" will form, which, in effect, will spontaneously break the local symme- 
try of the theory. Thus, the electric flux tube will become the boundary of a domain 
wall. As usual, this domain wall will decay by quantum tunneling — an electric flux 
tube will spontaneously nucleate, and expand, consuming the wall. Thus, there will 
be no stable electric flux tubes, and no inflnite range Aharonov-Bohm interaction 
between flux tubes and monopoles. 

We can check whether this expectation is correct in strong-coupling perturbation 
theory. We proceed by expanding e~'^p'''^ in powers of (3 at each plaquette, and e~^''"^'' 
in powers of A at each cube. Roughly speaking, the terms that survive when the C/;'s 

and rjp^s are summed arc ones such that a set of "tiled" cubes forms a closed (three- 
dimensional) hypersurface, or else the tiled cubes form an open hypersurface that is 
bounded by a (two-dimensional) surface of tiled plaquettes. In other words, strong- 
coupling perturbation theory can be interpreted as a sum over histories for (heavily 
suppressed) domain walls bounded by electric flux tubes. 

Consider, now, the behavior of G*^-^)(E), where (/) denotes the deflning represen- 
tation of SU{N). If /3 = 0, then the leading contribution to ^G'(-/")(S)^, for A « 1, 
is obtained by tiling the minimal hypersurface that is bounded by S; thus we have 

^G(-^)(E)^ - (^)Volume(S) _ ^^Q^^^g^, 

The interpretation is that G^-^\T,) is the boundary of a domain wall, where the wall 
tension is cr ~ — ln(l/A) (e is the lattice spacing). But for > 0, this domain 
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wall is unstable. When the surface S is very large, a much larger contribution to 
is obtained by tiling the plaquettes of S; this contribution is 

^^(■^^(E)^ - {p/N)^'''<^^ . (10.17) 

(As /? — > 0, the electric flux tube becomes arbitrarily heavy, and the domain wall is 
arbitrarily long-lived.) 

Accordingly, the operator W^^^C,!^) does not create a stable electric flux tube. 
The leading behavior of its expectation value, too, is found by tiling the plaquettes 
of E, so that 

(W^-^\C, E)^ - (/3/Ar)Area(S) _ (IQ.IS) 

There is no dependence on the area of the minimal surface bounded by C (compare 
eq. (7.2)), signifying that the dynamical flux tube created by W^^\C, E) has vanishing 
tension. 

It is obvious that the leading contribution to {^G(-^)(E)^ is unaffected by an 
insertion of Bn{C), even if C and E link. So we have 

lim (^eI!\e,C)^^1; (10.19) 

there is no long-range Aharonov-Bohm interaction. 

Incidentally, we could have chosen the cube action to be 

^ - X{r^^)^ + c.c . (10.20) 

Then, in effect, for A << 1, charge-m monopolcs condense; this breaks the local 
symmetry to Zm, where M is the greatest common factor of N and m (compare 
Section 4). In other words, W^'^\C,T,) does create a stable flux tube, if a source 
transforming as the representation (i/) has no Aharonov-Bohm interaction with a 
charge-m monopole, and if {u) represents Zm nontrivially. 
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10.5. Magnetic Hair 



Now we consider the regime (3 « 1 and A >> 1. In this hmit, virtual monopoles 
are heavily suppressed. We anticipate that quantum-mechanical magnetic charge 
fluctuations will not wipe out the inflnite-range Aharonov-Bohm interaction between 
monopoles and electric flux tubes. 

Again, we can check this expectation against perturbation theory. Weak-coupling 
perturbation theory, for A >> 1 is an expansion in the number of frustrated cubes. 
The frustrated cubes form closed loops that we may interpret as the world lines of 
virtual monopoles. The expansion in /3, as before, is an expansion in the number 
of tiled plaquettes. The tiled plaquettes form closed surfaces that we may interpret 
as the world sheets of electric flux tubes. Perturbation theory, then, is a sum over 
histories for (heavily suppressed) magnetic monopoles and electric flux tubes. 

For example, consider (^G^-^^T,)^ . The leading nontrivial contribution arises from 
a conflguration such that (in a particular gauge) rjp — e^'^'^^^^ for a single plaquette 
P contained in E*, while 77P = 1 for all other plaquettes. Flipping one plaquette 
frustrates four cubes, so this contribution gives 

(g(^)(E))~1+ J2 (e^-^/^-l) 

n=±l 

Summing disconnected contributions causes the result to exponentiate; we flnd 

^^(•^^(E)^ ~exp(-«('^^'^)Area(S)) , (10.22) 

where 

e^K^'^""^ - 2 (1 - cos(27r/iV)) exp (^-8X 

(and e is the lattice spacing). The interpretation is that, because of the Aharonov- 
Bohm interaction between monopoles and flux tubes, inserting G'^-'^-'(S) modifies the 
contribution to the vacuum energy due to virtual monopole pairs that wind around 
E, resulting in a renormalization of the tension of the classical string source. 



-2A 1 - cos 



2n 



(10.21) 



l-cos( - 



;i0.23) 
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When we compute (^W^-^\C, a similar renormalization of the tension of the 
classical string on E occurs. But in addition, for /5 << 1, the configurations that 
contribute have a surface of tiled plaquettes bounded by C. The leading behavior, 
then, is 

(w^f\C,E)^ ~ (/3/iV)^'^^^(^) exp (-K('^^'^)Area(S)) . (10.24) 
We conclude that W^^\C, E) creates a stable electric flux tube with string tension 

e2^(dyn) ^ -ln{/3/N) . (10.25) 

When the operator Bn{C) is inserted, it tends to frustrate the cubes in C*. But, 
as we already noted in Section 10.3, frustrated cubes can be avoided if the rjp^s assume 
a suitable configuration. We may choose an arbitrary surface T on the dual lattice 
whose boundary is C. Dual to the plaquettes of T is a set of plaquettes T* of the 
original lattice. The desired configuration (in a particular gauge) is 

rjp = e2^*^/^ , P e T* , 

(10.26) 

77P= 1 , P^T*. 

(The local symmetry transformation eq. (10.3) deforms the surface T, but leaves its 
boundary intact.) 

By summing over gauge field fiuctuations about the configuration eq. (10.26), we 
find the leading behavior 

(Pn(C)) - exp (^-Mi''^''^ Perimeter(C)) , (10.27) 

where 

eMt""^ r. 2N\P/Nf (^1 - cos ) (10.28) 

is the renormalization of the mass of the classical monopole source. This renormal- 
ization is associated with virtual electric flux tubes whose world sheets link the world 
line of the monopole, and arises because of the Aharonov-Bohm interaction between 
monopole and fiux tube. 
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When the operator Bn{Pij) is inserted, a hne of frustrated cubes connecting i 
and j cannot be avoided, and so we find the leading behavior 

{Bn{Pi,j)) - exp (-Mi'^^") Distance(i, j)) exp (-M^°") Length(P)) , (10.29) 

where 

eMi^^'^)~2A(^l-cos(^^^^ . (10.30) 

Evidently, for A » 1 and /3 « 1, we have Mi^^""^ » M^^''^ Indeed, (S„(C)) is 
dominated by small fluctuations about the configuration eq. (10.26) for precisely this 
reason — the rcnormalization of the classical source is much less costly than screening 
the source with dynamical monopoles. 

With the 77p's in the configuration eq. (10.26), the operator G(-^)(E) assumes the 
value exp{27rink / N) , where k is the linking number of E and C. Furthermore, except 
on the loop C, this configuration is locally equivalent to the trivial configuration with 
r]p — 1 everywhere. Thus, as we expand in the small fluctuations about eq. (10.26), 
we find, to each order of the expansion, 

hm (Ei/^) = [e^^^-Z^y^'"'''^ . (10.31) 

We see that, at least to each order of perturbation theory, our expectation is con- 
firmed. In a confining theory that contains weakly coupled dynamical magnetic 
monopoles, there is an infinite-range Aharonov-Bohm interaction between monopoles 
and electric fiux tubes. 
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APPENDIX 



In this appendix we consider in more detail some of the lattice perturbation 
theory calculations mentioned in the body of the text. Specifically, we calculate, in a 
pure gauge theory, the behavior of {Fa(T,,xo)) and (-Ba(S, C,xo)) in leading order in 
both strong and weak coupling perturbation theory. We demonstrate that (-Fa(S, xq)) 
exhibits an area law decay in both limits, and that only in the weak coupling limit 
does {Ba{T,,C,xo)) create a dynamical string on C. We also consider the problem 
that arises when one attempts to use the untraced Wilson loop operator eq. (5.15) to 
construct an Aharonov-Bohm order parameter. 

A.l. Pure gauge theory: strong coupling 

We start with the calculation of {Fa{T,,xo)) and C, xq)) in the strong 

coupling regime {P « 1) of the pure gauge theory defined by the plaquette action 
eq. (5.11), namely 

42ge,P = -/^x(^^(^P)+C-C- (A.l) 

For definiteness we will assume that the representation (R) that defines the theory 
is irreducible, and that it satisfies the constraints that R <Si R does not contain the 
trivial representation, while R<^ R* contains it exactly once. 

The expectation value (Fa(E,xo)) is given by 

(Fa(E,xo)) = ij^^ , (A.2) 

where the unnormalized expectation value of an operator X is defined by 

m) = 11 f (^n«^p(-'^^))- 

To find the renormalized string tension we need to calculate the leading behavior 
of both ((Fa(E, Xq))) and ((1)). The perturbation expansion in the strong coupling 
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regime is of the form of a sum over closed surfaces, in which surfaces of greater area are 
suppressed by powers of /? compared to smaller ones. Formally, the strong coupling 
expansion proceeds by performing a character expansion on the exponentiated Wilson 
action. 

The link integrations then select out closed surfaces S, each formed by "tiling" S with 
factors of x^'^\Up) for each plaquette P & S. The leading non-trivial contribution 
will be from the smallest surfaces, which furthermore are tiled with characters whose 
associated factors of C^^\(5) have the lowest non-trivial dependence on (3. The factors 
C('')(/3) are found by multiplying eq. (A.4) by x''''*\Up), summing U p over the group, 
and using the character orthogonality relations along with the assumed properties of 
the representation (i?), giving 

C^R){(3) = 1 + /3 + 0{f), C^^*\(3) = 1 + /3 + 0((3^), (A.5) 

with all other irreducible representations acquiring their first non-trivial /3 dependence 
at 0{P^). Therefore the leading contribution will come from surfaces tiled with 
characters in the representations (R) and (R*). 

In the evaluation of (Fa) we must, in addition, identify the leading contribution 
that is sensitive to the presence of F,,,. and hence does not cancel between {{Fa}) and 
((!)). All closed surfaces that do not intersect S are unaffected by the presence of 
Fa- Therefore we factorize both {{Fa}) and ((1)) in eq. (A. 2) into contributions from 
surfaces "on- E" and "off-E", the off-E dependence cancelling between them. The 
dominant on-E contribution is shown in Fig. 26, where all six faces are tiled with 
factors of either x^^^ or x^^ ^- This results in a contribution to {{Fa)) from Fig. 26of 

2((3/nR)'\x^^\a)\\ (A.6) 

where ur is the dimension of the representation (R). The factor of {nji)~^ comes 
from the six non-trivial link integrations, and the factor of 2 from the two equal 
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contributions of either x^-^) or \ Since the number of possible positions of this 
elementary cube on the surface S is given by the area of S (measured in terms of 
the size of the elementary plaquette e^), to get the contribution of a single elementary 
cube located anywhere on E we should further multiply (A. 6) by A{Y1). 

The area law decay of (F^) arises from considering an arbitrary number of such 
elementary tiled cubes on the surface E, rather than just a single contribution. In 
the "dilute gas" approximation where we can ignore correlations between the various 

elementary cubes, this leads to an exponentiation of the above result with a overall 
factor of A{Tj). Finally, the "on S" contributions to ((1)) lead to a similar exponential 
dependence on A(S), except that the group element a is replaced by the identity 
element. Putting these together yields the result 

(F«(E,a;o)) =exp(2(/3/n^)'(|x(^)(a)|2-n|)A(S)) . (A.7) 

In other words, the renormalized classical string tension defined in eq. (6.4) is given 

by 

to lowest order in /3. 

The calculation of (i?a(S, C, xq)) proceeds in a very similar way; the only question 
concerns the behavior at the boundary curve C. Here there exist contributions to 
(Sa(E, C, Xq)) similar to those of Fig. 26except only one plaquette P e E* is contained 
in the cube. This leads to an additional decay of (Sa(E, C, xq)) depending on the 
length P{C) of the perimeter (eq. (7.3)), 

{Ba{^, C, xo)) ~ exp A(E)) exp (-m^")p(C)) , (A.9) 

with 

em^r''^ (4 - nRRex^^\a)) . (A.IO) 

Since the configurations that dominate the expectation values of Fa and Ba do 
not extend deep into the volume enclosed by E, we find that there are no stable 



(A.8) 
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dynamical strings (or vortices) in this regime, and therefore that there is no long- 
distance Aharonov-Bohm interaction between flux tubes and charges. 

A. 2. Weak coupling 

This is the regime in which we expect to find A^a\^,xo;C) displaying a non- 
trivial Aharonov-Bohm interaction. To start we need to discover which configurations 
dominate the behavior of (Fo(E, xq)). When /3 >> 1, frustrated plaquettes (ones 
with non-minimal gauge action) are highly suppressed. However, as described in 
Section 6.2, an insertion of Fa{T,,xo) tends to frustrate the plaquettes in E*, unless 
the link variables assume a suitable configuration. The configuration that leaves 
no frustrated plaquettes, and thus acts as the "ground state" in the presence of Fa 
is illustrated in Fig. 9. It consists of a "forest" of links inside S that have value 
Ui = a. Specifically we choose a three-dimensional hypersurface Q made up of cubes 
of the dual lattice that has boundary E. Dual to these cubes is a set of links fl* on 
the original lattice. The configuration with no excited plaquettes is then given by 
eq. (6.7), namely Ui — a ior I e fl* , and Ui — e otherwise. The area law decay of 
(Fa(E, xq)) in the weak coupling regime arises from fiuctuations of the link variables 
around this state, which are sensitive to the presence of the non-trivial background. 

The leading such contribution to the decay of (Fa(E,a;o)) in the pure gauge the- 
ory is shown in Fig. 27a. (For simplicity we discuss the situation for vortices in three 
Euclidean dimensions.) It consists of a configuration that has nine excited plaquettes, 
and is constructed by considering two neighboring plaquettes in E*. These plaquettes 
are connected with four links with Ui — b to make a cube, with b summed over the 
group. In this configuration, there are eight excited plaquettes with plaquette action 
proportional to f3Kex{b) (the unshaded plaquettes in Fig. 27b). But, in addition, 
there is an excited plaquette (if a and b do not commute) with action proportional 
to pilexidba^^b^^) (the shaded plaquette in Fig. 27b). This extra excitation distin- 
guishes the contributions of Fig. 27 to ((Fa(E, xq))) and ((1)). Changes of variable 
can move the ninth excited plaquette around, but cannot get rid of it. 
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Of course, the process illustrated in Fig. 27 is only one of a set of similar processes 
where instead of summing over the same value b for each of the four connecting links, 
we sum independently over their values bi. These differ from the contribution of 
Fig. 27 in two ways: there will either be a greater number of excited plaquettes leading 
to a subdominant contribution, or, if some of the hi are taken to be the indentity, 
changes of variable on the remaining, independent, hi show that the a-links have no 
physical effect. This is why, for instance, we do not consider the less suppressed, and 
apparently non- trivial, contribution arising from the excitation of just two joining 
links. 

We can understand this in physical terms. The contribution to {Fa(T,, xq)) reflects 
how the non-trivial boundary condition introduced by the vortex (or string) affects 
the quantum fluctuations of the gauge fields. To see the effect, we must consider 
processes in which virtual gauge field excitations (with gauge quantum numbers that 
Aharonov-Bohm scatter off the inserted string of flux a) propagate around the string. 
The process illustrated in Fig. 27 is the leading one of this type. 

The end result of these considerations is a renormalization of the classical string 
tension (or vortex mass in three dimensions) leading to an area law (respectively, 
perimeter law) decay of (Fa(E, xq)) as in eq. (6.4). An estimate of the renormalized 
vortex mass in the weak-coupling regime arising from Fig. 27 is 



where we require more information about the group G to explicitly evaluate the sum. 

Also, note that wc have found area law decay of (Fa(S,a;o)) in both regimes of 
the pure non-abclian gauge theory described by cq. (A.l) (when a is not in the center 
of G). This differs from the situation in a pure abclian gauge theory where a change 
of variables shows that {Fa{T,,xo)) = ll^^ This is to be expected since the pure non- 



(ren) 




(A.ll) 



(l - exp{{-2(3)[nR - Rex^'^Hbab-'a-')]}) 
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abelian gauge theory contains excitations that, via the Aharonov-Bohm effect, are 
sensitive to the introduction of flux on E. 

To see that there exist stable dynamical strings in this regime of the pure gauge 
theory, we must consider (Sa(E, C, xq)), rather than xq)). Since the surface E 

now has a boundary, the three dimensional hypersurface Q must also end on a two 
dimensional surface S. This surface is dual to a set of plaquettes S* on the original 
lattice that must be frustrated, since only one of their links is contained within Q*. 
Clearly the dominant configuration is the one in which S is the minimal area surface 
with boundary C. Each such plaqucttc is suppressed relative to the corresponding 
vacuum configuration that dominates ((1)) by an amount 

exp{2p{Rex^'^\a)-nR)). (A.12) 

In other words (i?o(E, C, xq)) behaves as in eq. (7.2) 

(5«(E,C,xo)) = exp(-«:^'^^)^(E)) exp{-K^^'''^A{C)), (A.13) 

where A{C) is the area of the minimal surface S, and is the dynamical string 

tension of the stable string with fiux in the a conjugacy class, 

e'K'^^'^^^2i3{nR-Rex^^\a)). (A.14) 

In some situations it is possible that Sa(E, C, xq) creates not just one dynamical 
string, of fiux a, on C, but instead creates a number of strings of varying fiuxes 
ai, . . . ,ak, where aia2 . . .a^ = a. The corresponding link configuration would be a 
"forest state" of links of value Ui = a that terminates gradually in a number of steps 
rather than dropping to the indentity across a single plaquette. This is the more 
favorable configuration if the total area law suppression of the multiple strings is less 
than that for the single string. This translates into the condition (in the pure gauge 
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case, and ignoring perimeter corrections that vanish in the infinite size hmit) 

k 

{k - l)nR + Rex^^\a) - ^ Rex^^\ai) < . (A.15) 

i=l 

As an example of this, consider the expectation value, in a pure Zq gauge theory, of 
the operator Bs{T,, C, xq) that attempts to create a dynamical string with flux a = 
on C (where cu — exp(27ri/6)). The dominant process in this case is not the creation 
of an string, but instead that of three uj strings, since, substituting into eq. (A.15), 
we find 2 + (-1) - 3(1/2) < 0. 

We now turn to the order parameter A''a\Ti,xo;C), which is defined (eq. (5.3)) 
in terms of Fa(T,) and the Wilson loop operator W^'^\C). Fa{T,) has been described 
above, so we will now discuss the Wilson loop operator, and demonstrate eq. (5.4). 

First, let us consider the untraced Wilson operator U^'^\C,xo) (eq. (5.15)). It 
is easy to calculate the leading behavior of the expectation value of U^'^\C,xq) at 
weak coupling: as above we factorize the contributions from links on-C and off-C, 
and cancel the off-C contributions, 

{{U(-\C,xo))) _ Et/r--Et7.-P^'^H^i---^L)niiexp(4L>/3Rex(^)(^0) 
^^^^^ " E[/r--Ec/,aiiexp(4Z^/3Rex(^)(t/z)) 

(A.16) 

where we work in D + 1 Euclidean dimensions, and f/j for i = 1 . . . L are the links in 
C. (In eq. (A.16), we sum over all values of the link variables contained in C, but 
keep all other link variables fixed at Ui = 1.) Because we have not traced the Wilson 
loop the multiple summations factorize. We then use 



-Lj2D('''\g)x^^'\g) = (A.17) 
where I^'^^ is the identity matrix. Along with eq. (A. 4) this imphes that 



(aM(C,.„)).(^i!il(2M)%H. (A,8) 
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If we repeat the calculation in the background created by a string insertion operator 
Fa(E, xq) then the only difference is that for each time C links S there is a link in C 
that is multiplied by a. Neither Fa(E,a;o) nor Xq) is gauge-invariant, but a 

gauge transformation conjugates both a and U^'^^ by the value of the transformation 
at Xq, so these conjugations cancel, and the gauge-invariant result is that I^'^^ is 
replaced by D^'^\a'^^^''^^), where A;(E, C) is the linking number of the two curves. 
This proves eq. (5.16), and by taking a trace it proves eq. (5.4). 

However there is a major caveat that limits the usefulness of the untraced Wilson 
loop: there are corrections to eq. (5.16) that do not go to zero as C and S become 
large and well separated. Recall that Fa(E,a;o) corresponds to a thought experiment 
in which a specific group element a (not just a conjugacy class) was associated with a 
string, using an arbitrary basis at a point xq. As described in Section 5.2 this means 
that in defining Fa(E,xo) we give each plaquette in E* a long "tail" of links that 
stretches back to xq. The problem manifests itself (in a particular gauge) when we 
take into account configurations in which some (but not all) of the links coming out 
of xq arc excited. In these configurations the inserted flux may be conjugated, but 
not the measured group element (or vice versa). 

Concretely, the leading contribution, in weak coupling perturbation theory, that 
causes problems for the untraced Wilson loop is simply the excitation of a single link 
Ui = b on the path that connects xq to E. In the three dimensional case this leads 
to the excitation of four plaquettes, and is thus suppressed relative to the ground 
state configuration by a factor of exp(8/?(Re x'-^H^) ~ ^ii))- However, the action of 
these configurations, though large at weak coupling, is completely independent of the 
size and separation of E and C, and so remains finite as E and C become large and 
far apart. Their effect on the untraced order parameter is to conjugate the inserted 
fiux (and thus the fiux measured by U^'^\C,xq)). Therefore, there are corrections to 
eq. (5.16) that occur at finite order in weak coupling perturbation theory that render 
the measured flux uncertain up to conjugation. 

The physical interpretation of these conflgurations is that the path connecting xq 
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to S is linked by a virtual vortex-antivortex pair of flux b. In four dimensions there 

exist analogous processes where a virtual cosmic string links the connecting path. 
These effects do not correct the leading behavior of the ABOP given by eq. (5.4), be- 
cause the traced Wilson loop is gauge-invariant, and so is not sensitive to conjugation 
of the inserted flux. 
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FIGURE CAPTIONS 

1) (a) Standard paths ai and a2, based at xq, that enclose vortex 1 and vortex 
2. (b) Vortex 1 winds around vortex 2. (c) Path that, during the winding of 
vortex 1 around vortex 2, gets dragged to ai. (d) Path that gets dragged to a2 
during the winding. 

2) A particle transforming as the representation (z/) of H interacts with a vortex 
pair. The particle and its antiparticle are initially prepared in an uncharged 
(gauge-singlet) state. After the particle winds around the vortex, the particle- 
antiparticle pair has acquired a nontrivial charge. 

3) (a) Standard paths a and f3, based at xq, associated with the a string and the 
b string, (b) Path that is dragged to the path (3 as the loop of b string winds 
through the loop of a string (in the sense defined by a) . 

4) The branching of a c = a6 string into an a string and a b string. The standard 
paths, based at xq, associated with each string are indicated. 

5) (a) The standard paths a, (3, and 7 associated with entangled strings. These 
paths are mapped by the string homomorphism to group elements a, b, and 
c. (b) The paths /? and are both mapped to b, if the bridge connecting the 
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strings can be removed by crossing the strings, (c) Standard paths associated 
with the "upper halves" of entangled strings, (d) The paths [3 and f3" are both 
mapped to b, if the bridge can be removed by crossing the strings "the other 
way." 

6) (a) Standard paths that are mapped by the string homomorphism to elements 
a, 6 e G that are not in the unbroken group H. Each string is the boundary of 
a domain wall, (b) A path that is mapped to ab^^ e H. This path does not 
cross the wall. 

7) (a) Two holes appear in a domain wall, due to the spontaneous nucleation of 
loops of a string and b string, (b) The holes meet and coalesce. A segment of 
ab~^ string appears that bridges the hole, (c) The b string has higher tension 

than the a string, and so pulls the string junctions toward each other, (d) A 
loop of ab~^ string breaks free from the decaying wall. 

8) A vein in a wall. 

9) The closed curve S in 2+1 dimensions. The dashed line is S, comprised of links 
of the dual lattice. The plaqucttcs shown are those in E*, which are dual to 
the links of E. The links marked by arrows are those in Q*; they are dual to 
the plaquettes in a surface Q that is bounded by E. (We may also interpret the 
dashed line as a shce through the closed surface E, in 3+1 dimensions.) 

10) A slice through a closed surface E that links once with the loop C. C intersects 
the hypersurface fl once. 

11) (a) Spacetime diagram (in 2+1 dimensions) showing a virtual b vortex- 
anti vortex pair that nucleates, winds around the basepoint xq, and reanni- 
hilates. If the path based at xq at time ti is assigned the flux a, then the path 
at time t2 is assigned flux bab^^. (b) The same process, but with xq now taken 
to be a fixed point in Euclidean spacetime. The flux of the classical a vortex 
is measured to be bab~^, due to the effect of the virtual b vortex pair. Shaded 
areas are surfaces bounded by the vortex world lines. 



79 



12) Paths Pi, P2, -P3 used in the construction of the operator 

-^01,02,03(^1; ^2, S3; Xq). 

13) A different choice for the path P2 from the basepoint xq to the world sheet E2. 

14) The paths Pi and P2, based at xq, and the open curves Ei and E2 used in 
the definition of the 't Hooft operator 5^ 5(2:1, yi, Ei, 0:2, 1/2, E2, xq), in 2+1 
dimensions. 

15) Two topologically inequivalent ways of choosing the hypersurface fli that is 
bounded by Ei. 

16) Surface of minimal area bounded by the loops Ci and C2 (when the loops are 
close together). 

17) The classical world sheet Edass and the dynamical world sheet E^yn associated 
with the 't Hooft loop operator Pa(C', ^class; ^^o)- 

18) Two contributions to the mixing of an a string and a b string. In (a), an ab~^ 
string spontaneously nucleates, expands, meets the a string, and converts it 
into a b string. In (b), the a string shrinks and annihilates, then the b string 
nucleates and grows. 

19) Linked world hues Ei and E2 of a classical a vortex and a classical b vortex. If 
a and b do not commute, then the classical vortices must exchange a dynamical 
aba~^b~^ vortex. The exchange occurs along the path P, the shortest path that 
connects the two world hues. (Compare Fig. 5.) 

20) A plaquette (shaded) contained in the intersection of fll and fl^. 

21) World sheets Ei and E2 of a classical a string and classical b string that intersect 
at two isolated points. (Only a shce through E2 is shown.) The indicated shaded 
region is the world sheet of a dynamical aba~^b~^ string that connects the two 
classical strings. 

22) World sheets Ei and E2 of a classical a string and a classical b string. (Only a 
slice through E2 is shown.) Ei is a torus that links with E2. The shaded region 
is the world sheet of a dynamical aba~^b~^ string. 
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23) The Borromean rings. Si is the world hne of an a vortex, E2 is the world line 

of a 6 vortex, and C is the world line of a charged particle that transforms as 
the representation (z/). The charged particle transfers charge to the a vortex- 
antivortex pair, and the charge is subsequently detected via the Aharonov- 
Bohm interaction of the pair with the b vortex. 

24) Nonintersecting surfaces Qi and Q2 that are bounded by the vortex world hues 
El and E2. 

25) The path P from y to x that is used to construct the gauge-invariant correlation 
function K^f^\x,y, P). A classical charged particle propagates along P, and a 
dynamical charged particle propagates from x to y. 

26) The leading contribution to (Fa(E, xq)), in the strong-coupling limit of the pure 
gauge theory, arises from tiling a cube (shaded) that contains two plaquettes of 
E*. 

27) The leading contribution to (Fa(E, xq)), in the weak-coupling hmit of the pure 
gauge theory. The cube shown contains two plaquettes that are in E*, and the 
link variables assume the indicated values. (Unmarked links have the value Ui — 
1.) There are altogether nine excited plaquettes — eight (unshaded) plaquettes 
with Up = b and one (shaded) plaquette with Up = aba~^b~^. 
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